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We propose a universal, on-chip quantum transducer based on surface acoustic waves in piezo¬ 
active materials. Because of the intrinsic piezoelectric (and/or magnetostrictive) properties of the 
material, our approach provides a universal platform capable of coherently linking a broad array of 
qubits, including quantum dots, trapped ions, nitrogen-vacancy centers or superconducting qubits. 

The quantized modes of surface acoustic waves lie in the gigahertz range, can be strongly confined 
close to the surface in phononic cavities and guided in acoustic waveguides. We show that this type 
of surface acoustic excitations can be utilized efficiently as a quantum bus, serving as an on-chip, 
mechanical cavity-QED equivalent of microwave photons and enabling long-range coupling of a wide 
range of qubits. 


I. INTRODUCTION 

The realization of long-range interactions between re¬ 
mote qubits is arguably one of the greatest challenges to¬ 
wards developing a scalable, solid-state spin based quan¬ 
tum information processor [1]. One approach to address 
this problem is to interface qubits with a common quan¬ 
tum bus which distributes quantum information between 
distant qubits. The transduction of quantum information 
between stationary and moving qubits is central to this 
approach. A particularly efficient implementation of such 
a quantum bus can be found in the field of circuit QED 
where spatially separated superconducting qubits inter¬ 
act via microwave photons confined in transmission line 
cavities [2-4]. In this way, multiple qubits have been cou¬ 
pled successfully over relatively large distances of the or¬ 
der of millimeters [5, 6]. Fueled by the dramatic advances 
in the fabrication and manipulation of nanomechanical 
systems [7], an alternate line of research has pursued the 
idea of coherent, long-range interactions between indi¬ 
vidual qubits mediated by mechanical resonators, with 
resonant phonons playing the role of cavity photons [8- 
13]. 

In this paper, we propose a new realization of a quan¬ 
tum transducer and data bus based on surface acous¬ 
tic waves (SAW). SAWs involve phonon-like excitations 
bound to the surface of a solid and are widely used in 
modern electronic devices e.g. as compact microwave fil¬ 
ters ]14, 15]. Inspired by two recent experiments ]16, 17], 
where the coherent quantum nature of surface acoustic 
waves (SAWs) has been explored, here we propose and 
analyze SAW phonon modes in piezo-active materials as 
a universal mediator for long-range couplings between 
remote qubits. Our approach involves qubits interacting 
with a localized SAW phonon mode, defined by a high¬ 
er resonator, which in turn can be coupled weakly to a 
SAW waveguide serving as a quantum bus; as demon¬ 
strated below, the qubits can be encoded in a great vari¬ 


ety of spin or charge degrees of freedom. We show that 
the Hamiltonian for an individual node (for a schematic 
representation see Fig.l) can take on the generic Jaynes- 
Gummings form {h = 1), 

■f^node = TWca'I'a-I-g (cr+a-I-cr“a'l’) , (1) 

where a refers to the usual Pauli matrices describing the 
qubit with transition frequency Wq and a is the bosonic 
operator for the localized SAW cavity mode of frequency 
ujc 127r GHz ]18]. The coupling g between the qubit and 
the acoustic cavity mode is mediated intrinsically by the 



Figure 1: (color online). SAW as a universal quantum trans¬ 
ducer. Distributed Bragg reflectors made of grooves form an 
acoustic cavity for surface acoustic waves. The resonant fre¬ 
quency of the cavity is determined by the pitch p, fc = Vsl2p. 
Reflection occurs effectively at some distance inside the grat¬ 
ing; the fictitious mirrors above the surface are not part of 
the actual experimental setup, but shown for illustrative pur¬ 
poses only. Red arrows indicate the relevant decay channels 
for the cavity mode: leakage through the mirrors, internal 
losses due to for example surface imperfections, and conver¬ 
sion into bulk modes. Qubits inside and outside of the solid 
can be coupled to the cavity mode. In more complex struc¬ 
tures, the elastic medium can consist of multiple layers on top 
of some substrate. 
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piezo-properties of the host material, it is proportional 
to the electric or magnetic zero-point fluctuations associ¬ 
ated with a single SAW phonon and, close to the surface, 
can reach values of 5 ~ 400MHz, much larger than the 
relevant decoherence processes and sufficiently large to 
allow for quantum effects and coherent coupling in the 
spin-cavity system as evidenced by cooperativities [19] of 
C ~ 10 — 100 [see Section IV and Tab. I for definition 
and applicable values]. For Wq « Uc, iFnode allows for 
a controlled mapping of the qubit state onto a coherent 
phonon superposition, which can then be mapped to an 
itinerant SAW phonon in a waveguide, opening up the 
possibility to implement on-chip many quantum commu¬ 
nication protocols well known in the context of optical 
quantum networks [13, 20] . 

The most pertinent features of our proposal can be 
summarized as follows: (1) Our scheme is not specific 
to any particular qubit realization, but—thanks to the 
plethora of physical properties associated with SAWs 
in piezo-active materials (strain, electric and magnetic 
fields)—provides a common on-chip platform accessible 
to various different implementations of qubits, compris¬ 
ing both natural (e.g., ions) and artificial candidates such 
as quantum dots or superconducting qubits. In partic¬ 
ular, this opens up the possibility to interconnect dis¬ 
similar systems in new electro-acoustic quantum devices. 
(2) Typical SAW frequencies lie in the gigahertz range, 
closely matching transition frequencies of artificial atoms 
and enabling ground state cooling by conventional cryo¬ 
genic techniques. (3) Our scheme is built upon an estab¬ 
lished technology [14, 15] : Lithographic fabrication tech¬ 
niques provide almost arbitrary geometries with high pre¬ 
cision as evidenced by a large range of SAW devices such 
as delay lines, bandpass Alters or resonators etc. In par¬ 
ticular, the essential building blocks needed to interface 
qubits with SAW phonons have already been fabricated, 
according to design principles familiar from electromag¬ 
netic devices: (i) SAW resonators, the mechanical equiva¬ 
lents of Fabry-Perot cavities, with low-temperature mea¬ 
surements reaching quality factors of Q ~ 10 ^ even at gi¬ 
gahertz frequencies [21-23], and (ii) acoustic waveguides 
as analogue to optical fibers [14] . (4) For a given fre¬ 
quency in the gigahertz range, due to the slow speed of 
sound of approximately ^ lO^m/s for typical materials, 
device dimensions are in micrometer range, which is con¬ 
venient for fabrication and integration with semiconduc¬ 
tor components, and about 10 ® times smaller than cor¬ 
responding electromagnetic resonators. (5) Since SAWs 
propagate elastically on the surface of a solid within a 
depth of approximately one wavelength, the mode vol¬ 
ume is intrinsically confined in the direction normal to 
the surface. Further surface confinement then yields large 
zero-point fluctuations. ( 6 ) Yet another inherent advan¬ 
tage of our system is the intrinsic nature of the coupling. 
In piezoelectric materials, the SAW is accompanied by an 
electrical potential (j) which has the same spatial and tem¬ 
poral periodicities as the mechanical displacement and 
provides an intrinsic qubit-phonon coupling mechanism. 


For example, recently qubit lifetimes in GaAs singlet- 
triplet qubits were found to be limited by the piezoelec¬ 
tric electron-phonon coupling [24] . Here, our scheme pro¬ 
vides a new paradigm, where coupling to phonons be¬ 
comes a highly valuable asset for coherent quantum con¬ 
trol rather than a liability. 

In what follows, we first review the most important 
features of surface acoustic waves, with a focus on the 
associated zero-point fluctuations. Next, we discuss the 
different components making up the SAW-based quan¬ 
tum transducer and the acoustic quantum network it 
enables: SAW cavities, including a detailed analysis of 
the achievable quality factor Q, SAW waveguides and a 
variety of different candidate systems serving as qubits. 
Lastly, as exemplary application, we show how to trans¬ 
fer quantum states between distant nodes of the network 
under realistic conditions. Finally, we draw conclusions 
and give an outlook on future directions of research. 

II. SAW PROPERTIES 

Elastic waves in piezoelectric solids are described by 

pili Cijjildjdl'U,}^ — ( 2 ) 

CijOidjcj) — Cijf^didf^Uj (3) 

where the vector u (x, t) denotes the displacement held 
(x is the cartesian coordinate vector), p is the mass den¬ 
sity and repeated indices are summed over (f, j = x, y, z); 
c, € and e refer to the elasticity, permittivity and piezo¬ 
electric tensors, respectively [25] ; they are largely defined 
by crystal symmetry [26]. For example, for cubic crys¬ 
tals such as GaAs there is only one non-zero component 
for the permittivity and the piezoelectric tensor, labeled 
as e and eu, respectively [25] . Since elastic disturbances 
propagate much slower than the speed of light, it is com¬ 
mon practice to apply the so-called quasi-static approxi¬ 
mation [26] where the electric field is given by = —di(j). 
When considering surface waves, Eq.(2) and (3) must 
be supplemented by the mechanical boundary condition 
that there should be no forces on the free surface (taken 
to be at z = 0 with z being the outward normal to the 
surface), that is T^x = T^y = T^z = 0 at z = 0 (where 
Tij = CijkidiUk + ekijdkfj) is the stress tensor), and appro¬ 
priate electrical boundary conditions [25] . 

If not stated otherwise, the term SAW refers to the 
prototypical (piezoelectric) Rayleigh wave solution as 
theoretically and experimentally studied for example in 
Refs. [16, 17, 25, 27] and used extensively in different elec¬ 
tronic devices [14, 15] . It is non-dispersive, decays ex¬ 
ponentially into the medium with a characteristic pene¬ 
tration depth of a wavelength and has a phase velocity 
Vs = uj/k that is lower than the bulk velocities in that 
medium, because the solid behaves less rigidly in the ab¬ 
sence of material above the surface [26]. As a result, it 
cannot phase-match to any bulk-wave [14, 28] . As usual, 
we consider specific orientations for which the piezoelec¬ 
tric held produced by the SAW is strongest [14, 28], for 
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example a SAW with wavevector along the [110] direction 
of a (001) GaAs crystal [cf. Refs. [16, 25[ and Appendices 
A and B[. 

SAWs in quantum regime. —In a semiclassical picture, 
an acoustic phonon associated with a SAW creates a 
time-dependent strain field, su = {diUk + dkUi)/2 and 
a (quasi-static) electrical potential (j) (x, t). Upon quan¬ 
tization, the mechanical displacement becomes an op¬ 
erator that can be expressed in terms of the elemen¬ 
tary normal modes as u(x) = [v„ (x) a„ -I- h.c.j, 

where a„ (aj^) are bosonic annihilation (creation) op¬ 
erators for the vibrational eigenmode n and the set 
of normal modes v„ (x) derives from the Helmholtz- 
like equation >Vv„ (x) = —pa;^v„ (x) associated with 
Eq.(2) and (3). The mode normalization is given by 
J" d^xpv* (x) • v„ (x) = h/2LOn [24, 29[. An important 
figure of merit in this context is the amplitude of the me¬ 
chanical zero-point motion Uq. Along the lines of cavity 
QED [2[, a simple estimate for Uq can be obtained by 
equating the classical energy of a SAW ~ f d^xpxi^ with 
the quantum energy of a single phonon, that is hu. This 
leads directly to 

Uo « ^hl2pVsA, (4) 

where we used the dispersion relation lo = Vgk and the 
intrinsic mode confinement V « AX characteristic for 
SAWs. The quantity Uq refers to the mechanical ampli¬ 
tude associated with a single SAW phonon close to the 
surface. It depends only on the material parameters p 
and Vs and follows a generic ~ A~^l'^ behaviour, where 
A is the effective mode area on the surface. The estimate 
given in Eq.(4) agrees very well with more detailed calcu¬ 
lations presented in Appendix C. Several other important 
quantities which are central for signal transduction be¬ 
tween qubits and SAWs follow directly from Uq'. The (di¬ 
mensionless) zero-point strain can be estimated as sq ~ 
kUo- The intrinsic piezoelectric potential associated with 
a single phonon derives from Eq.(3) as (po « (ei 4 /e) Uq 
[30[. Lastly, the electric field amplitude due to a single 
acoustic phonon is ~ k(j)Q = (en/e) kUo, illustrating 
the linear relation between electric field and strain char¬ 
acteristic for piezoelectric materials [8[. In summary, we 
typically find Uq ~ 2fm/ y/A [/rm^], yielding Uq « 2fm for 
micron-scale confinement (cf. Appendix D) . This is com¬ 
parable to typical zero-point fluctuation amplitudes of lo¬ 
calized mechanical oscillators [31[. Moreover, for micron- 
scale surface confinement and GaAs material parameters, 
we obtain ~ 20V/m which compares favorably with 
typical values of ~ 10“^V/m and ~ 0.2V/m encountered 
in cavity and circuit QED, respectively [2[. 

For the sake of clarity, we have focused on piezo- 
electric materials so far. However, there are also piezo- 
magnetic materials that exhibit a large magnetostrictive 
effect. In that case, elastic distortions are coupled to a 
(quasi-static) magnetic instead of electric field [32, 33[; 
for details see Appendix D. For typical materials such 
as Terfenol-D the magnetic field associated with a single 
phonon can be estimated as Hg « (2 — 6)p,T/-\/A 



Figure 2: (color online). Characterization of a groove-based 
SAW cavity, (a) Quality factor Q for A = 100 (dashed blue) 
and N = 300 (red solid) grooves as a function of the normal¬ 
ized grove depth h/Xc- For shallow grooves, Q is limited by 
leakage losses due to imperfect acoustic mirrors (Qr-iegime, 
gray area), whereas for deep grooves conversion to bulk modes 
dominates (Qi,-regime); compare asymptotics (dash-dotted 
lines), (b) Ratio of desired to undesired decay rates Kgd/Kbd- 
The stronger Q is dominated by Qr, the higher fCgd/^bd. Here, 
w/p = 0.5, D = 5.25Ac, and /c = 3GHz; typical material pa¬ 
rameters for LiNbOa have been used [cf. Appendix E[. 

Finally, we note that composite structures comprising 
both piezoelectric and piezomagnetic materials can sup¬ 
port magneto-electric surface acoustic waves [34, 35[. 

III. SAW CAVITIES & WAVEGUIDES 

SAW cavities .—To boost single phonon effects, it is es¬ 
sential to increase Ug. In analogy to cavity QED, this can 
be achieved by confining the SAW mode in an acoustic 
resonator. The physics of SAW cavities has been theoret¬ 
ically studied and experimentally verified since the early 
1970s [14, 36[. Here, we provide an analysis of a SAW 
cavity based on an on-chip distributed Bragg reflector 
in view of potential applications in quantum information 
science; for details see Appendix E. SAW resonators of 
this type can usually be designed to host a single res¬ 
onance fc = Wc/27r = Vs/Xc (Ac = 2p) only and can 
be viewed as an acoustic Fabry-Perot resonator with ef¬ 
fective reflection centers, sketched by localized mirrors 
in Fig.l, situated at some effective penetration distance 
into the grating [14[. Therefore, the total effective cav¬ 
ity size along the mirror axis is L^. > D, where D is 
the physical gap between the gratings. The total cavity 
line-width n = Wc/Q = ^gd + ^bd can be decomposed 
into desired (leakage through the mirrors) and undesired 
(conversion into bulk modes and internal losses due to 
surface imperfections etc.) losses, labeled as Kgd and 
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Kbd, respectively; for a schematic illustration compare 
Fig.l. For the total quality factor Q, we can typically 
identify three distinct regimes [cf. Fig.2]: For very small 
groove depths h/Xc < 2%, losses are dominated by cou¬ 
pling to SAW modes outside of the cavity, dubbed as 
Qr-regime (Kgd 3> Kbd), whereas for very deep grooves 
losses due to conversion into bulk-modes become exces¬ 
sive (Qb-regime, Kgd ^ ^bd)- In between, for a suf¬ 
ficiently high number of grooves N, the quality factor 
Q can ultimately be limited by internal losses (surface 
cracks etc.), referred to as Qm-regime (Kgd ^ ^bd)- For 
N « 300, we find that the onset of the bulk-wave limit 
occurs for h/Xc > 2.5%, in excellent agreement with ex¬ 
perimental findings [36, 37]. With regard to applica¬ 
tions in quantum information schemes, the Qj.-regime 
plays a special role in that resonator phonons leaking 
out through the acoustic mirrors can be processed fur¬ 
ther by guiding them in acoustic SAW waveguides (see 
below). To capture this behaviour quantitatively, we 
analyze /Cgd/^bd [cf. Fig.2[: for Kgd/nhd ^ !> leakage 
through the mirrors is the strongest decay mechanism for 
the cavity phonon, whereas the undesired decay channels 
are suppressed. Our analysis shows that, for gigahertz 
frequencies /c « 3GHz, N « 100 and h/Xc ~ 2%, a 
quality factor of Q ~ 10^ is achievable, together with an 
effective cavity confinement Lc « 40Ac (for D < 5Ac) 
and Kgd/Kbd ^ 20 (illustrated by the circle in Fig.2); 
accordingly, the probability for a cavity phonon to leak 
through the mirrors (rather than into the bulk for exam¬ 
ple) is Kgd/ (^gd + Kbd) ^ 95%. Note that the resulting 
total cavity linewidth of K/27r = fc/Q ~ (1 — 3) MHz is 
similar to the ones typically encountered in circuit QED 
[6[. To compare this to the effective cavity-qubit cou¬ 
pling, we need to fix the effective mode area of the SAW 
cavity. In addition to the longitudinal confinement by the 
Bragg mirror (as discussed above) a transverse confine¬ 
ment length Ltrans (in direction y) can be provided, e.g., 
using waveguiding, etching or (similar to cavity QED) 
focusing techniques [14, 38, 39[. For transverse confine¬ 
ment Ttrans ~ (1 ~ 5) /xm and a typical resonant cavity 
wavelength Ac ~ l^m, the effective mode area is then 
A = itransAc ~ (40 — 200)/im^. In the desired regime 
Kgd/Khd ^ 1) this is largely limited by the deliberately 
low reflectivity of a single groove; accordingly, the cavity 
mode leaks strongly into the grating such that Lc ^ D 
[cf. Appendix E for details]. While we have focused on 
this standard Bragg design (due to its experimentally 
validated frequency selectivity and quality factors), let 
us shortly mention potential approaches to reduce A and 
thus increase single-phonon effects even further: (i) The 
most straightforward strategy (that is still compatible 
with the Bragg mirror design) is to reduce Ac as much 
as possible, down to the maximum frequency fc = Vs/Xc 
that can still be made resonant with the (typically highly 
tunable) qubit’s transition frequency ujql2TT', note that 
fundamental Rayleigh modes with fc « 6GHz have been 
demonstrated experimentally [40] . (ii) In order to in¬ 
crease the reflectivity of a single groove, one could use 


deeper grooves. To circumvent the resulting increased 
losses into the bulk [cf. Fig. 2(b)], free-standing struc¬ 
tures (where the effect of bulk phonon modes is reduced) 
could be employed, (iii) Lastly, alternative cavity designs 
such as so-called trapped energy resonators make it pos¬ 
sible to strongly confine acoustic resonances in the center 
of plate resonators [41] . 

SAW waveguides .—Not only can SAWs be confined in 
cavities, but they can also be guided in acoustic waveg¬ 
uides (WGs) [14, 42] . Two dominant types of design are: 

(i) Topographic WGs such as ridge-type WGs where the 
substrate is locally deformed using etching techniques, or 

(ii) overlay WGs (such as strip- or slot-type WGs) where 
one or two strips of one material are deposited on the 
substrate of another to form core and clad regions with 
different acoustic velocities. If the SAW velocity is slower 
(higher) in the film than in the substrate, the film acts as 
a core (cladding) for the guide whereas the unmodified 
substrate corresponds to the cladding (core). An atten¬ 
uation coefficient of ~ 0.6dB/mm has been reported for 
a lO^m-wide slot-type WG, defined by A1 cladding lay¬ 
ers on a GaAs substrate [38, 39] . This shows that SAWs 
can propagate basically dissipation-free over chip-scale 
distances exceeding several millimeters. Typically, one¬ 
dimensional WG designs have been investigated, but— 
to expand the design flexibility—one could use multiple 
acoustic lenses in order guide SAWs around a bend [28] . 


IV. UNIVERSAL COUPLING 

Versatility .—To complete the analogy with cavity 
QED, a non-linear element similar to an atom needs to 
be introduced. In the following, we highlight three dif¬ 
ferent exemplary systems, illustrating the versatility of 
our SAW-based platform. We focus on quantum dots, 
trapped ions and NV-centers, but similar considerations 
naturally apply to other promising quantum information 
candidates such as superconducting qubits [7, 8, 17, 43], 
Rydberg atoms [44] or electron spins bound to a phos¬ 
phorus donor atom in silicon [40] . In all cases consid¬ 
ered, a single cavity mode a, with frequency Wc close to 
the relevant transition frequency, is retained. We provide 
estimates for the single-phonon coupling strength and co- 
operativity [cf. Tab.I], while more detailed analyses go 
beyond the scope of this work and are subject to future 
research. 

(i) QD Gharge qubit: A natural candidate for our 
scheme is a charge qubit embedded in a lithographically 
defined GaAs double quantum dot (DQD) containing a 
single electron. The DQD can well be described by an 
effective two-level system, characterized by an energy off¬ 
set e and interdot tunneling tc yielding a level splitting 
n = V Lt1 [45] . The electron’s charge e couples to 
the piezoelectric potential; the deformation coupling is 
much smaller than the piezoelectric coupling and can 
therefore safely be neglected [46] . Since the quantum 
dot is small compared to the SAW wavelength, we ne- 
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charge qubit (DQD) 

spin qubit (DQD) 

trapped ion 

NV-center 

coupling g 

(200 - 450)MHz 

(10 - 22.4)MHz 

(1.8-4.0)kHz 

(45 - 101)kHz 

cooperativity C 

11 - 55 

21 - 106 

7-36 

10-54 


Table I: Estimates for single-phonon coupling strength g and cooperativity C. We have used ^4 = (1 — 5)fj,m x 40Ac, T = 20mK 
[17], (conservative) quality factors of Q = (1,1,3,1) x 10® and frequencies of ujc = 27r(6,1.5,2x 10 ®,3) GHz for the four 
systems listed. For the spin qubit = 2/rs [68], and for the trapped ion scenario, gion(C'ion) is given for d = 150/rm due to the 
prolonged dephasing time further away from the surface (Cion improves with increasing d, even though gion decreases, up to a 
point where other dephasing start to dominate). Further details can be found in the main text. 


gleet potential effects coming from the structure making 
up the dots (heterostructure and metallic gates); for a de¬ 
tailed discussion see Appendix I. Performing a standard 
rotating-wave approximation (valid for <5, gch ^ ^c)i we 
find that the system can be described by a Hamiltonian 
of Jaynes-Cummings form, 

Of 

HAot = {S^a +S~a^) , (5) 

where 5 = — Wc specifies the detuning between the 

qubit and the cavity mode, and = |±) (=f| (and 
so on) refer to pseudo-spin operators associated with 
the eigenstates |±) of the DQD Hamiltonian (cf. Ap¬ 
pendix F). The Hamiltonian Ffdot describes the coher¬ 
ent exchange of excitations between the qubit and the 
acoustic cavity mode. The strength of this interaction 
(/ch = ecjjoT {kd)sui{kl/2) is proportional to the charge 
e and the piezoelectric potential associated with a sin¬ 
gle phonon (pQ. The decay of the SAW mode into the 
bulk is captured by the function (kd) [d is the dis¬ 
tance between the DQD and the surface; see Appendix 
B for details], while the factor sm{kl/2) reflects the as¬ 
sumed mode function along the axis connecting the two 
dots, separated by a distance 1. For (typical) values 
of I « Ac/2 = 250nm and d « 50nm <C Ac, the geo¬ 
metrical factor F {kd) sm-{kl/2) then leads to a reduc¬ 
tion in coupling strength compared to the bare value 
600 (at the surface) by a factor of ^ 2 only. In total, 
we then obtain ^ch ~ 2GHz/Ay^4~[/rnW[. For lateral con¬ 
finement Ttrans ~ (1 ~ 5) /im, the effective mode area is 
^ = -^transAc ~ (20 — 100)/rm^. The resulting charge- 
resonator coupling strength ^ch ~ (200 — 450)MHz com¬ 
pares well with values obtained using superconducting 
qubits coupled to localized nano-mechanical resonators 
made of piezoelectric material where g « (0.4 — 1.2)GHz 
[7, 8] or superconducting resonators coupled to Gooper 
pair box qubits ( 5/271 « 6MHz) [3], transmon qubits 
( 5/277 « lOOMHz) [47] and indium arsenide DQD qubits 
( 5/277 « 30MHz) [48]. Note that, in principle, the cou¬ 
pling strength 5ch could be further enhanced by addition¬ 
ally depositing a strongly piezoelectric material such as 
LiNbOa on the GaAs substrate [16]. Moreover, with a 
LiNbOa film on top of the surface, also non-piezoelectric 
materials such as Si or Ge could be used to host the 
quantum dots [49]. The level splitting fl (t) and interdot 
tunneling tc {t) can be tuned in-situ via external gate 
voltages. By controlling S one can rapidly turn on and 
off the interaction between the qubit and the cavity: For 


an effective interaction time r = 77/2539 (geff = 2gchtc/^) 
on resonance ((5 = 0), an arbitrary state of the qubit is 
swapped to the absence or presence of a cavity phonon, 
i.e., (a|-)+/3|-k))|0) ^ |-) (a |0) - 70 |1)), where |n) 
labels the Fock states of the cavity mode. Apart from this 
SWAP operation, further quantum control techniques 
known from cavity QED may be accessible [50]. Note 
that below we will generalize our results to spin qubits 
embedded in DQDs. 

(ii) Trapped ion: The electric field associated with the 
SAW mode does not only extend into the solid, but, for 
a free surface, in general there will also be an electrical 
potential decaying exponentially into the vacuum above 
the surface ~ exp[—fc|z|] [25]; cf. Appendix B. This 
allows for coupling to systems situated above the surface, 
without any mechanical contact. For example, consider a 
single ion of charge q and mass m trapped at a distance 
d above the surface of a strongly piezoelectric material 
such as LiNbOa or AIN. The electric dipole induced by 
the ion motion couples to the electric field of the SAW 
phonon mode. The dynamics of this system are described 
by the Hamiltonian 

Hion = ujea^a + LUtb^b -I- 5ion {ab^ + a%) , (6) 

where b refers to the annihilation operator of the ion’s 
motional mode and ujt is the (axial) trapping frequency. 
The single phonon coupling strength is given by 5ion = 
qxo X kepoF (ked) = qpoF (ked) Apart from the 

exponential decay F {kd) = exp [—kd], the effective cou¬ 
pling is reduced by the Lamb-Dicke parameter 5 ld = 
277 x 0 /Ac, with xq = •\/hl2muJt, since the the motion of 
the ion is restricted to a region small compared with 
the SAW wavelength Ac- For LiNbOa, a surface mode 
area of A = (1 — 5)/im x 40Ac, the commonly used ^Be'*’ 
ion and typical ion trap parameters with d « 30/im and 
074/277 « 2MHz [51], we obtain 5ion ~ (3 —6.7)kHz. Here, 
5ion refers to the coupling between the ion’s motion and 
the cavity. However, based on iJion, one can in princi¬ 
ple generalize the well-known protocols operating on the 
ion’s spin and motion to operations on the spin and the 
acoustic phonon mode [52]. 

(iii) NV-center: Yet another system well suited for our 
scheme are NV centers in diamond. Even though dia¬ 
mond itself is not piezoactive, it has played a key role 
in the context of high-frequency SAW devices due to 
its record-high sound velocity [14]; for example, high- 
performance SAW resonators with a quality factor of 
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Figure 3: (color online), (a) Schematic illustration for cou¬ 
pling to a NV center via a piezo-magnetic (PM) material (see 
text for details); surface grooves (not shown) can be used to 
provide SAW phonon confinement, (b) SAWs can be gen¬ 
erated electrically based on standard interdigital transducers 
(IDTs) deposited on the surface. Typically, an IDT consists 
of two thin-film electrodes on a piezo-electric material, each 
formed by interdigitated fingers. When an ac voltage is ap¬ 
plied to the IDT on resonance (defined by the periodicity of 
the fingers as ojidt/Stt = Vs/piot, where Va is the SAW prop¬ 
agation speed), it launches a SAW across the snbstrate surface 
in the two directions perpendicular to IDT fingers[14, 15, 17]. 

Q = 12500 at Wc ^ lOGHz were experimentally demon¬ 
strated for AlN/diamond heterostructures [53, 54]. To 
make use of the large magnetic coupling coefficient of the 
NV center spin = Stt x 28GHz/T, here we consider a 
hybrid device composed of a thin layer of diamond with 
a single (negatively charged) NV center with ground- 
state spin S implanted a distance d « lOnm away from 
the interface with a strongly piezo-magnetic material. 
Equivalently, building upon current quantum sensing ap¬ 
proaches [55, 56[, one could use a diamond nanocrystal 
(typically ~ lOnm in size) in order to get the NV centre 
extremely close to the surface of the piezo-magnetic ma¬ 
terial and thus maximize the coupling to the SAW cavity 
mode; compare Fig. 3(a) for a schematic illustration. In 
the presence of an external magnetic field Bext [57[, the 
system is described by 

Hj^v = DS^ + 7NvBext • S -I- LOcCt^a (7) 

+ffNV J 

cic—x,y,z 

where D = 27rx2.88GHz is the zero-field splitting, (^nv = 
7nvBo is the single phonon coupling strength and is 
a dimensionless factor encoding the orientation of the NV 
spin with respect to the magnetic stray field of the cavity 
mode. For d <C Ac, a rough estimate shows that at least 
one component of ry^y is of order unity [33[. For a NV 
center close to a Terfenol-D layer of thickness h ^ Xc, 
we find ^nv ~ 400kHz /A [/rm^]. Thus, as compared to 
direct strain coupling < 200Hz /A [pm^], the presence 
of the piezomagnetic layer is found to boost the single 
phonon coupling strength by three orders of magnitude; 


this is in agreement with previous theoretical results for 
a static setting [33[. 

Decoherence. —In the analysis above, we have ignored 
the presence of decoherence which in any realistic set¬ 
ting will degrade the effects of coherent qubit-phonon 
interactions. In this context, the cooperativity param¬ 
eter, defined as C = g^T 2 / [^(hth + 1)], is a key figure 
of merit. Here, T 2 refers to the corresponding dephasing 
time, while nth = (exp [huJc/ksT] — 1)”^ gives the ther¬ 
mal occupation number of the cavity mode at tempera¬ 
ture T. The parameter C compares the coherent single¬ 
phonon coupling strength g with the geometric mean of 
the qubit’s decoherence rate ~ and the cavity’s ef¬ 
fective linewidth ~ k (nth + 1); in direct analogy to cav¬ 
ity QED, C > 1 marks the onset of coherent quantum 
effects in a coupled spin-oscillator system, even in the 
presence of noise; cf. Ref.[10[ and Appendix H for a de¬ 
tailed discussion. To estimate C, we take the following 
parameters for the dephasing time T 2 : For system (i) 
T 2 « 10ns has been measured close to the charge degen¬ 
eracy point e = 0 [45[. In scenario (ii) motional decoher¬ 
ence rates of 0.5Hz have been measured in a cryogenically 
cooled trap for an ion height of 150/im and IMHz mo¬ 
tional frequency [51[. Since this rate scales as ~ 

[52, 58[, we take T 2 [s] « 2(d[/xm] /150)^. Lastly, for the 
NV-center (iii) T 2 « 0.6s has been demonstrated for en¬ 
sembles of NV spins [59[ and we assume an optimistic 
value of T 2 = 100ms, similarly to Ref.[60[. The results 
are summarized in Tab.I. We find that C > 1 should be 
experimentally feasible which is sufficient to perform a 
quantum gate between two spins mediated by a thermal 
mechanical mode [9[. 

Qubit-qubit coupling. —When placing a pair of qubits 
into the same cavity, the regime of large single spin co¬ 
operativity C 1 allows for coherent cavity-phonon- 
mediated interactions and quantum gates between the 
two spins via the effective interaction Hamiltonian Hint = 
gdr{S^S 2 + h.c.), where g^r = g^/S ^ g in the so-called 
dispersive regime [4[. For the estimates given in Tab.I, we 
have restricted ourselves to the Q^-regime with Q « 10^, 
where leakage through the acoustic mirrors dominates 
over undesired (non-scalable) phonon losses (Kgd S> Kbd)- 
However, note that small-scale experiments using a sin¬ 
gle cavity only (where there is no need for guiding the 
SAW phonon into a waveguide for further quantum in¬ 
formation processing) can be operated in the Q^-regime 
(which is limited only by internal material losses), where 
the quality factor Q « Qm ^ 10® is maximized (and thus 
overall phonon losses minimal). 

As a specific example, consider two NV-centers, both 
coupled with strength ^nv ~ lOOkHz to the cavity and 
in resonance with each other, but detuned from the res¬ 
onator. Since for large detuning 6 the cavity is only virtu¬ 
ally populated, the cavity decay rate is reduced to Kdr = 
{g^K « « IkHz (for /c = 3GHz, Q = 2 x 10®), 

whereas the spin-spin coupling is g^v ~ 0.1(;nv ~ lOkHz. 
Therefore, T 2 = 1ms is already sufficient to approach the 
strong-coupling regime where gdr 7^ Kdr,72~^- 
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Finally, we note that, in all cases considered above, one 
could implement a coherent, electrical control by pump¬ 
ing the cavity mode using standard interdigital transduc¬ 
ers (IDTs) [14, 15, 17]; compare Fig. 3(b) for a schematic 
illustration. The effect of the additional Hamiltonian 
Hdrive = 5 cos (widtO [a + Can be accounted for by 
replacing the cavity state by a coherent state, that is 
a —>■ a. For example, in the case of Eq.(5), one could 
then drive Rabi oscillations between the states |-|-) and 
|—) with the amplified Rabi frequency fl/j = ga. 


V. STATE TRANSFER PROTOCOL 

The possibility of quantum transduction between 
SAWs and different realizations of stationary qubits en¬ 
ables a variety of applications including quantum infor¬ 
mation achitectures that use SAW phonons as a quan¬ 
tum bus to couple dissimilar and/or spatially separated 
qubits. The most fundamental task in such a quantum 
network is the implementation of a state transfer protocol 
between two remote qubits 1 and 2, which achieves the 
mapping (a |0)^ -f/3 |1) J(g)[ 0)2 -t |0)^(g)(a [ 0)2 + /3 11 ) 2 )- 
In analogy to optical networks, this can be accomplished 
via coherent emission and reabsorption of a single phonon 
in a waveguide [13]. As first shown in the context of 
atomic QED [20], in principle perfect, deterministic state 
transfer can be implemented by identifying appropriate 
time-dependent control pulses. 

Before we discuss a specific implementation of such a 
transfer scheme in detail, we provide a general approxi¬ 
mate result for the state transfer fidelity As demon¬ 
strated in detail in Appendix H, for small infidelities one 
can take 

T ^l-e-CC-\ ( 8 ) 

as a general estimate for the state transfer fidelity. Here, 
individual errors arise from intrinsic phonon losses ~ 
£ = Khd/Kgd and qubit dephasing ~ C~^ ~ re¬ 

spectively; the numerical coefficient C ^ O (1) depends 
on the specific control pulse and may be optimized for 
a given set of experimental parameters [61] . This sim¬ 
ple, analytical result holds for a Markovian noise model 
where qubit dephasing is described by a standard pure 
dephasing term leading to an exponential loss of coher¬ 
ence ^ exp(—t/T 2 ) and agrees well with numerical re¬ 
sults presented in Ref.[61[. For non-Markovian qubit de¬ 
phasing an even better scaling with C can be expected 
[9[. Using experimentally achievable parameters e « 5% 
and C « 30, we can then estimate T « 90%, showing 
that fidelities sufficiently high for quantum communica¬ 
tion should be feasible for all physical implementations 
listed in Tab. I. 

In the following, we detail the implementation of a 
transfer scheme based on spin qubits implemented in 
gate-defined double quantum dots (DQDs) [62] . In par¬ 
ticular, we consider singlet-triplet-like qubits encoded in 



Figure 4: (color online). (a) Average fidelity IF of 
the state transfer protocol for a coherent superposition 
IV') = (|0) “ |1)) /\/2 in the presence of quasi-static (non- 
Markovian) Overhauser noise, as a function of the root- 
mean-square fluctuations Cnuc in the detuning parameters 
Si {i = 1,2), for fCbd/Kgd = 0 (solid line, circles) and 

Kbd/^gd = 10% (dash-dotted line, squares), (b) After n = 100 
runs with random values for Si, T approximately reaches con¬ 
vergence. The curves refer to (Tnuc/^gd = (0, 2,..., 10)% 
(from top to bottom) for Kbd/w^gd = 10%. (c) Pulse shape 
g\ (t) for first node. 


lateral QDs, where two electrons are localized in adja¬ 
cent, tunnel-coupled dots. As compared to the charge 
qubits discussed above, this system is known to feature 
superior coherence timescales [63-66]; these are largely 
limited by the relatively strong hyperfine interaction be¬ 
tween the electronic spin and the nuclei in the host en¬ 
vironment [64], resulting in a random, slowly evolving 
magnetic (Overhauser) field for the electronic spin. To 
mitigate this decoherence mechanism, two common ap¬ 
proaches are (i) spin-echo techniques which allow to ex¬ 
tend spin-coherence from a time-ensemble-averaged de¬ 
phasing time T 2 « 100ns to T 2 > 250/rs [66], and (ii) 
narrowing of the nuclear field distribution [64, 67] . Re¬ 
cently, real-time adaptive control and estimation meth¬ 
ods (that are compatible with arbitrary qubit operations) 
have allowed to narrow the nuclear spin distribution to 
values that prolong to > 2/rs [68] . For our pur¬ 
poses, the latter is particularly attractive as it can be 
done simply before loading and transmitting the quan¬ 
tum information, whereas spin-echo techniques can be 
employed as well, however at the expense of more com¬ 
plex pulse sequences (see Appendix for details). In or¬ 
der to couple the electric field associated with the SAW 
cavity mode to the electron spin states of such a DQD, 
the essential idea is to make use of an effective elec¬ 
tric dipole moment associated with the exchange-coupled 
spin states of the DQD [69-72]. As detailed in Appendix 















G, we then find that in the usual singlet-triplet sub¬ 
space spanned by the two-electron states {I'flJJ-) , IJJ-tl)}, 
a single node can well be described by the prototypical 
Jaynes-Cummings Hamiltonian given in Eq.(l). As com¬ 
pared to the direct charge coupling ^chj the single phonon 
coupling strength g is reduced since the qubit states \l) 
have a small admixture of the localized singlet (<S'o 2 lO 
{I = 0,1) only. Using typical parameters values, we find 
g « 0.1(?ch ~ 200MHz/-y A [gTo?] [73]. In this system, the 
coupling g (t) can be tuned with great flexibility via both 
the tunnel-coupling tc and/or the detuning parameter e. 

The state transfer between two such singlet-triplet 
qubits connected by a SAW waveguide can be adequately 
described within the theoretical framework of cascaded 
quantum systems, as outlined in detail for example in 
Refs. [13, 20, 74, 75[: The underlying quantum Langevin 
equations describing the system can be converted into an 
effective, cascaded Master equation for the system’s den¬ 
sity matrix p. For the relevant case of two qubits, it can 
be written as p = £ideaiP + UnoiseP) where 


-^idealP — ^ Hg (t) -\- 7/tgd , p 

-|-2Kgd27 [oi -|- 02 ] p, (9) 

>C„oiseP = 2«;bd ^ 'D [Ui] p-i'^Si [S'/, p] . (10) 


2 = 1,2 


Here, I? [a] p = apa^ — | {a^a,p} is a Lindblad term 
with jump operator a and Hg (t) = Hi (t), with 
Hi{t) = gi{t)[S^ai -I- S'~a|] describes the coherent 
Jaynes-Cummings dynamics of the two nodes. The ideal 
cascaded interaction is captured by Uideai which contains 
the non-local coherent environment-mediated coupling 
transferring excitations from qubit 1 to qubit 2 [76[, while 
'Cnoise Summarizes undesired decoherence processes: We 
account for intrinsic phonon losses (bulk-mode conversion 
etc.) with a rate Kbd and (non-exponential) qubit de¬ 
phasing. Since the nuclear spins evolve on relatively long 
time-scales, the electronic spins in quantum dots typi¬ 
cally experience non-Markovian noise leading to a non¬ 
exponential loss of coherence on a characteristic time- 
scale T 2 given by the width of the nuclear field dis¬ 
tribution (Tnuc as T 2 = V^/fTnuc [M, 68 [. Recently, a 
record-low value of CT„uc/27r = 80kHz has been reported 
[ 68 [, yielding an extended time-ensemble-averaged elec¬ 
tron dephasing time of T 2 = 2.8/iS. In our model, to 
realistically account for the dephasing induced by the 
quasi-static, yet unknown Overhauser field, the detuning 
parameters Si are sampled independently from a normal 
distribution p (Si) with zero mean (since nominal reso¬ 
nance can be achieved via the electronic control param¬ 
eters) and standard deviation dnuc [67[; see Appendix 
G for details. In Appendix J we also provide numer¬ 
ical results for standard Markovian dephasing, showing 
that non-Markovian noise is beneficial in terms of faithful 
state transfer. 

Under ideal conditions where /Inoise = 0, the setup 
is analogous to the one studied in Ref.[20[ and the same 


time-symmetry arguments can be employed to determine 
the optimal control pulses gi (t) for faithful state trans¬ 
fer: if a phonon is emitted by the first node, then, upon 
reversing the direction of time, one would observe per¬ 
fect reabsorption. By engineering the emitted phonon 
wavepacket such that it is invariant under time reversal 
and using a time-reversed control pulse for the second 
node 52 (t) = gi (—t), the absorption process in the sec¬ 
ond node is a time-reversed copy of the emission in the 
first and therefore in principle perfect. Based on this 
reasoning (for details see Ref.[20[), we find the explicit, 
optimal control pulse shown in Fig. 4(c). 

To account for noise, we simulate the full master equa¬ 
tion numerically. The results are displayed in Fig. 4(a), 
where for every random pair S = (Si,S2) the fidelity of 
the protocol is defined as the overlap between the target 
state |'0tar) and the actual state after the transfer p(tf), 
that is Ts = (V’tar|p(^/) I'/’tar)- The average fidelity ^ of 
the protocol is determined by averaging over the classical 
noise in S, that is IF = f dSidS2P (Si) p (S2) SFi- Taking an 
effective mode area A k, 100/rm^ as above and Q « 10^ 
to be well within the Q^-regime where Kbd/^gd ~ 5%, 
we have g « Kgd ~ 20MHz. For two nodes separated 
by millimeter distances, propagation losses are negligi¬ 
ble and Atbd/wgd ~ 5% captures well all intrinsic phonon 
losses during the transfer. We then find that for realistic 
undesired phonon losses Kbd/^gd ~ 5% and crnuc/27r = 
80kHz (such that crnuc/%d ~ 2.5%) [68[, transfer fideli¬ 
ties close to 95% seem feasible. Notably, this could be 
improved even further using spin-echo techniques such 
that T 2 « 10^T2 [66[. Therefore, state transfer fideli¬ 
ties SF > 2/3 as required for quantum communication 
[77[ seem feasible with present technology. Near unit fi¬ 
delities might be approached from further optimizations 
of the system’s parameters, the cavity design, the control 
pulses and/or from communication protocols that correct 
for errors such as phonon losses [78-80] . Once the trans¬ 
fer is complete, the system can be tuned adiabatically 
into a storage regime which immunizes the qubit against 
electronic noise and dominant errors from hyperfine in¬ 
teraction with ambient nuclear spins can be mitigated by 
standard, occasional refocusing of the spins [66, 69] . Al¬ 
ternatively, one could also investigate silicon dots: while 
this setup requires a more sophisticated hetero-structure 
including some piezo-electric layer (as studied experimen¬ 
tally in Ref.[40[), it potentially benefits from prolonged 
dephasing times T 2 > 100p,s [81], since nuclear spins are 
largely absent in isotopically purified ^®Si. 


VI. SUMMARY & OUTLOOK 

In summary, we have proposed and analyzed SAW 
phonons in piezo-active materials (such as GaAs) as 
a universal quantum transducer that allows to convert 
quantum information between stationary and propagat¬ 
ing realizations. We have shown that a sound-based 
quantum information architecture based on SAW cavi- 
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ties and waveguides is very versatile, bears striking simi¬ 
larities to cavity QED and can serve as a scalable medi¬ 
ator of long-range spin-spin interactions between a vari¬ 
ety of qubit implementations, allowing for faithful quan¬ 
tum state transfer between remote qubits with existing 
experimental technology. The proposed combination of 
techniques and concepts known from quantum optics and 
quantum information, in conjunction with the technolog¬ 
ical expertise for SAW devices, is likely to lead to further, 
rapid theoretical and experimental progress. 

Finally, we highlight possible directions of research go¬ 
ing beyond our present work: First, since our scheme 
is not specific to any particular qubit realization, novel 
hybrid systems could be developed by embedding dissim¬ 
ilar systems such as quantum dots and superconducting 
qubits into a common SAW architecture. Second, our 
setup could also be used as a transducer between differ¬ 
ent propagating quantum systems such as phonons and 
photons. Light can be coupled into the SAW circuit via 
(for example) NV-centers or self-assembled quantum dots 
and structures guiding both photons and SAW phonons 
have already been fabricated experimentally [38, 39]. Fi¬ 
nally, the SAW architecture opens up a novel, on-chip 
test-bed for investigations reminiscent of quantum op¬ 
tics, bringing the highly developed toolbox of quantum 
optics and cavity-QFD to the widely anticipated field of 
quantum acoustics [11, 16, 17, 82[. Potential applica¬ 
tions include quantum simulation of many-body dynam¬ 
ics [83[, quantum state engineering (yielding for exam¬ 
ple squeezed states of sound), quantum-enhanced sens¬ 
ing, sound detection, and sound-based material analysis. 
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Appendix A: Classical Description of 
Nonpiezoelectric Surface Acoustic Waves 

In this Appendix, we review the general (classical) the¬ 
oretical framework describing SAW in cubic lattices, such 
as diamond or GaAs. We derive an analytical solution 
for propagation in the [110[ direction. The latter is of 
particular interest in piezoelectric systems. The classical 
description of SAW is explicitly shown here to make our 



Cll 

C 12 

C 44 

p[kg/m®] 

ei4 

Diamond 

107.9 

12.4 

57.8 

3515 

0 

GaAs 

12.26 

5.71 

6.00 

5307 

0.157 


Table II: Material properties [25[ for both diamond and GaAs. 
The elastic tensor c has three independent parameters, given 
in units of [10^°N/m^], while the piezoelectric tensor e has a 
single independent parameter ei 4 for cubic materials (units of 
C/mQ. 


work self-contained, but follows standard references such 
as Refs.[25[. 

Wave equation. —The propagation of acoustic waves 
(bulk and surface waves) in a solid is described by the 
equation 


ry I ■ ■ 

pu,(x,t) = |^, (Al) 

where u denotes the displacement vector with Ui being 
the displacement along the cartesian coordinate Xi {xi = 
x,X 2 = y,X 3 = z), p gives the mass density and T is 
the stress tensor; is the ith component of force per 
unit area perpendicular to the ij-axis. Moreover, x is 
the cartesian coordinate vector, where in the following 
we assume a material with infinite dimensions in i, y, 
and a surface perpendicular to the i-direction at z = 0 . 
The stress tensor obeys a generalized Hooke’s law (stress 
is linearly proportional to strain) 




(A2) 


where the strain tensor is defined as 


Ukl 


1 / ^ dui \ 

2 \ dxi dxk ) 


(A3) 


Using the symmetry Cijki = Cijik, in terms of displace¬ 
ments we find 


rp _ ^ 

Uj - Cyfci , 


(A4) 


such that Eq.(Al) takes the form of a set of three coupled 
wave equations 


pUj(x,t) Cijkl 


d'^Uk 

dxjdxi 


= 0 . 


(A5) 


The elasticity tensor c obeys the symmetries Cijki = 
Cjiki = c.,jik = Ckuj and is largely defined by the crys¬ 
tal symmetry. 

Mechanical boundary condition. —The free surface at 
z = 0 is stress free (no external forces are acting upon 
it), such that the three components of stress across z = 0 
shall vanish, that is T 13 = T 23 = r 33 = 0. This results in 
the boundary conditions 

Tiz = = 0 atz = 0. (A 6 ) 

0 x 1 
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k [7r/yum] 


Figure 5: (color online). Dispersion relation tOn = v„k of the 
three (n = 1, 2,3) Rayleigh-type SAW modes for propagation 
along x'W [110]. If not stated otherwise, we refer to the the 
lowest frequency solution as the SAW mode (solid line). 


(and cyclic permutations) while the mechanical boundary 
conditions can be written as 


^13 

723 

733 


C44 

C44 


/ du^ \ 

\dx ^ dz ) 

f dUy \ 


= 0 , 
= 0 , 


dU;, 

Cll ^-h C12 

az 


/ duy \ 

\ clx dy ) 


= 0 , 


(A8) 

(A9) 

(AlO) 


at 2 ; = 0. In the following we seek for solutions which 
propagate along the surface with a wavevector k = 
k {lx + my), where I = cos (0), m = sin (0) and 0 is the 
angle between the x-axis and k. Following Ref. [84], we 
make the ansatz 


Cubic lattice. —For a cubic lattice (such as GaAs or 
diamond) the elastic tensor Cijki has three independent 
elastic constants, generally denoted by cii,ci2, and C44; 
compare Tab. II. Taking the three direct two-fold axes as 
the coordinate axes, the wave equations then read 


( 1 = ( v\ (All) 

UJ Uj 


d'^Ux 


Cll 


d'^Ux 

dx'^ 


+ C 44 


d'^Ux d'^Ux 
dy'^ dz'^ 


+ (Cl 2 + C44) 


d'^Uy 

dxdy 


dxdz 


(A7) 


where the decay constant q describes the exponential de¬ 
cay of the surface wave into the bulk and c is the phase 
velocity. Plugging this ansatz into the mechanical wave 
equations can be rewritten as A4A = 0, where 


Clip + C 44 — 9^) — pc^ Im (ci2 + C44) 

M = \ Im (ci2 + C44) Ciim^ -|- C44 (P — q^) — pc^ 

Iq (ci2 + C44) mq (ci2 + C44) ciiq'^ - C44 -|- pc^ 


Iq (ci2 + C44) 

mq (ci2 + C44) 
..2 


(A12) 


and A = {U,V,iW). Nontrivial solutions for this homo¬ 
geneous set of equations can be found if the determinant 
of A4 vanishes, resulting in the so-called secular equation 
det (A4) = 0. The secular equation is of sixth order in q; 
as all coefficients in the secular equation are real, there 
are, in general, three complex-conjugate roots 
with the phase velocity c and propagation direction 0 as 
parameters. If the medium lies in the half space z > 0, 
the roots with negative real part will lead to a solution 
which does not converge as z —)■ 00 . Thus, only the roots 
which lead to vanishing displacements deep in the bulk 
are kept. Then, the most general solution can be written 
as a superposition of surface waves with allowed q^ values 


as 


{ux,uy,iu,)= Y. 

r^l,2,3 

(A13) 

where, for any q^ = q^ (c, 0), the ratios of the amplitudes 
can be calculated according to 


Ur _Vr _ iWr 
rjr Cr 


(AI4) 


where we have introduced the quantities 




Ciim^ -I- C44 (P — qr) — pc^ mqr (ci2 + C44) 
mqr (ci2 + C44) Ciig2 - C44 -h pc^ 


Vr 


mqr (ci2 + C44) Im (ci2 + C44) 
Cliqr - Cii + pc^ lqr{Cl2+Cii) 


(AI5) 
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and 


c. 


Im (ci2 + C44) Ciim^ + C44 [P — q ^) — pcP 

kr (ci2 + C44) mqr (ci2 + C44) 


(A 16 ) 

Note that for each root qr and displacement Ui there 
is an associated amplitude. The phase velocity c, how¬ 


ever, is the same for every root qr, and needs to be 
determined from the mechanical boundary conditions 
as described below. Similarly to the acoustic wave 
equations, the boundary conditions can be rewritten as 
B {Ki, K2, K3) = 0, where the boundary condition ma¬ 
trix B is 


B = 


( Kl — 9lCl K 2 — 92^2 Ks — 93^3 

m-Ci - 9i9i rnC2 - 92??2 wCs - 93% 

kl + mT]i + aqiCi ^^2 + W 7?2 + 092(2 ^(3+^-773-1-093(3 


(A 17 ) 


with a = Cii/ci 2 . Again, nontrivial solutions are found 
for det (S) = 0. The requirements det(Al) = 0, 
det (B) = 0 together with Eq.(A14) constitute the formal 
solution of the problem [84]; det (A4) = 0 and det (B) = 0 
may be seen as determining and q^, and Eq.(A14) then 
gives the the ratios of the components of the displace¬ 
ment. In the following, we discuss a special case where 
one can eliminate the 9 -dependence in det (B) — 0, lead¬ 
ing to an explicit, analytically simple equation for the 
phase velocity c, which depends only on the material 
properties. 

Propagation in [110] direction. —The wave equations 
simplify for propagation in high-symmetry directions. 
Here, we consider propagation in the [110]-direction, for 
which I = m = 1/V2; we define the diagonal as x' = 
{x + y) /\/2. Subtracting the second row from the first 
in A4, one finds that the common factor (cn — C 12 ) /2 — 
C 449 ^ — pc^ divides through the first row, which then be¬ 
comes (1,—1,0). Thus, U = V and the wave equations 
can be simplified to Mno {U,iW) = 0, where 


Miiii={ ~ ^ 

y ^29 (ci2 + C44) 0119^ - C44 -h pc^ 

(A18) 

with kl = (cii -|- C12 + 2C44) / 2 . Then, the secular equa¬ 
tion det (All 10) = 0 is found to be 

(c'li - pk - C449^) (c44 - pc^ - Cll9^) 

+ (ci 2 + C44) 9 ^ = 0,(A19) 



yielding the roots 94,9I. We choose the roots commensu¬ 
rate with the convergence condition yielding the general 
ansatz 



(A 20 ) 


with Ux = Uy = Ux' /v^. The amplitude ratios 7' = 
iWr/U'r can be obtained from the kernel of A 4 as 


ir = 9r 


C12 + C44 

C44 - Cii {X + q]) ’ 


(A 21 ) 


where X = pt? jcn. In the coordinate system {x',z}, the 
mechanical boundary conditions read 


duz 

dP 


C12 


dux' 

dx' 


dux 


Cll 


dz 

duz 


dz 


0, (z = 0) 

0. (z = 0) 


(A 22 ) 

(A 23 ) 


For the ansatz given in Eq.(A 20 ), they can be reformu¬ 
lated as Hiio {U[, U2) = 0 with 


^110 = 


li - 9i 72 - 92 I 

1 + 1 + ) ■ 


(A 24 ) 


The requirement det (-8410) = 0 can be written as 

9l [ci2 + pc ^ + Cll 9 i] [ci2 (c44 - pc ^) + C11C4492] “ 

92 [ci2 + pc^ + C1192] [ci2 (C44 “ PC^) + CiiC449l] = 0. 

From the symmetry of this equation it is clear that one 
can remove a factor (94 — 92) leading to 




+ C449i92 + C 12 X 



- C 449192 



0. 


Using simple expressions for 949I and q] + 9I obtained 
from Eq.(A 19 ), one arrives at the following explicit equa¬ 
tion for the wave velocity c [ 25 , 84 [ 



(A 25 ) 

which is cubic in A = pc^/cn. If not stated otherwise, we 
consider the mode with the lowest sound velocity, referred 
to as Rayleigh mode; compare Fig. 5 . 

Using the secular equation given in Eq.(A 19 ) and the 
mechanical boundary conditions, the ansatz given in 
Eq.(A 20 ) can be simplified as follows: The roots compat¬ 
ible with convergence in the bulk are complex conjugate. 
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H-/U 


Figure 6: (color online). Depth dependence of the (normal¬ 
ized) vertical displacement UzjU along x'\ \ [110] for a Rayleigh 
surface acoustic wave propagating on a (001) GaAs crystal. 
The acoustic amplitude decays away from the surface into the 
bulk on a characteristic length scale approximately given by 
the SAW wavelength X — 2n/k fa 1/rm. 


i.e. q = qi = q 2 , and therefore 7 = 7^ =72- Then, using 
the first row in the boundary condition matrix [compare 
Eq.(A24)], we can deduce 

U[ = C/e"*^, (A26) 

where 



7-9 


In summary, we find the following solution [25] 

u^, = [/(e-'?'=^-*‘^-hh.c.)e*'=(“'"=*), 

iu, = [/(7e-«'="-*‘^ + h.c.)e*'=("'-"‘), (A28) 

where the material-dependent parameters c, q, 7 and ip 
are determined by Eq.(A25), Eq.(A19), Eq.(A21) and 
Eq.(A25), respectively. For the GaAs parameters given 
in Tab. II, we get c = 2878m/s, g = 0.5 -I- 0.48z, 7 = 
—0.68+1.16i and ip = 1.05, respectively. The correspond¬ 
ing (normalized) transversal displacement is displayed in 
Fig. 6. 


Appendix B: Surface Acoustic Waves In 
Piezoelectric Materials 

In a piezoelectric material, elastic and electromagnetic 
waves are coupled. In principle, the field distribution can 
be found only by solving simultaneously the equations 
of both Newton and Maxwell. The corresponding solu¬ 
tions are hybrid elasto-electromagnetic waves, i.e., elas¬ 
tic waves with velocity Vg accompanied by electric fields, 
and electromagnetic waves with velocity c « lO^u^ ac¬ 
companied by mechanical strains. For the first type of 
wave, the magnetic field is negligible, because it is due to 


an electric field traveling with a velocity Vg much slower 
than the speed of light c; therefore, one can approxi¬ 
mate Maxwell’s equations as V x E = —d'B/dt « 0, 
giving E = — V(f>. Thus, the propagation of elastic waves 
in a piezoelectric material can be described within the 
quasi-static approximation, where the electric field is es¬ 
sentially static compared to electromagnetic fields [26[. 
The potential (f) and the associated electric field are not 
electromagnetic in nature but rather a component of the 
predominantly mechanical wave propagating with veloc¬ 
ity Vg. 


1. General Analysis 


Wave equation .—The basic equations that govern the 
propagation of acoustic waves in a piezoelectrical mate¬ 
rial connect the mechanical stress T and the electrical 
displacement D with the mechanical strain and the elec¬ 
trical field. The coupled constitutive equations are 


T ■ — 

-^7.1 - 


Dz = 


C^ijkl 


duk 

dxi 

d(j) 

'dxj 


^kij 


^ijk 


d(j) 

dxk 

duj 

dxk 


(Bl) 


where e with {cijk = eikj) and e are the piezoelectric and 
permittivity tensor, respectively. Here, Hooke’s law is 
extended by the additional stress term due to the piezo¬ 
electric effect, while the equation for the displacement Di 
includes the polarization produced by the strain. There¬ 
fore, Newton’s law becomes 


P'^i — ^ijkl 


d'^Uk 

dxjdxi 


“1“ ^kij Tx • 

UXjUXk 


(B2) 


For an insulating solid, the electric displacement Di must 
satisfy Poisson’s equation dDijdxi = 0 which yields 


d'^Uj d‘^(t) 

OXiOXk OXiOXj 

A(j) = 0, 


z> 0 
z > 0. 


(B3) 

(B4) 


Mechanical boundary conditions .—In the presence of 
piezoelectric coupling the mechanical boundary condi¬ 
tions [compare Eq.(A6)[ generalize to 

= Cizki^^ + ek^z-p^ =0 at z = 0. (B5) 

0X1 OXk 

Using the symmetries c^Jkl = Cjm and Ckij = Ckji it is 
easy to check that this is equivalent to Eq.(41) in Ref.[25[. 

Electric boundary condition .—In addition to the stress- 
free boundary conditions, piezoelectricity introduces an 
electric boundary condition: The normal component of 
the electric displacement needs to be continuous across 
the surface [38[, that is 

Dz {z = 0 +) =Dz{z = 0 -) , 


(B6) 
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where by definition = Czjkduj/dxk—ezjd(l)/dxj. Out¬ 
side of the medium (z < 0 ), we assume vacuum; thus, 
Dz = cqEz = —Codipout/dz, where the electrical poten¬ 
tial has to satisfy Poisson’s equation At/fout = 0. The 
ansatz 


(B7) 

gives A(/)out = (— (j)out = 0. Thus, for proper 
convergence far away from the surface z — )■ —oo, we take 
the decay constant 0 = 1 ; accordingly, the electrical po¬ 
tential decays exponentially into the vacuum above the 
surface on a typical length scale given by the SAW wave¬ 
length A = 27r/k ~ 1/rm. Therefore, for the electrical dis¬ 
placement outside of the medium, we find Dz = —eokcj). 
Lastly, the electrical potential has to be continuous across 
the surface [25], i.e., 

(j)(z = 0+) = (/fout (z = O') , (B 8 ) 

which allows us to determine the amplitude Aout- In 
summary, Eq.(B 6 ) ca be rewritten as 

{czjkduj/dxk - ezjd4>ldxj + eok(t>)\^^^ = 0. (B9) 



Figure 7: (color online). The dimensionless function J-{z) 
determines the decay of the electrical potential away from 
the surface into the bulk; the characterisitc length scale is 
approximately set by the SAW wavelength A = 27r/fe « 
Ifim. Inset: Density plot of the (normalized) electric poten¬ 
tial Re [<j)\ /4 )q = —T (kz) sin {kx' — uit) along x'\\ [110] for a 
Rayleigh surface acoustic wave propagating on a (001) GaAs 
crystal at t = 0. 


2. Cubic lattice 


Cubic lattice. —For a cubic, piezoelectric system there 
is only one independent nonzero component of the piezo¬ 
electric tensor called ei 4 [25, 26[. With this piezoelectric 
coupling, the wave equations are given by four coupled 
partial differential equations 


d'^Ux 


eAcj) 


Cll 


dx"^ 


+ C 44 


d'^Ux 


d'^Ux 


+ (Cl 2 + C44) 


d^Uy 

dxdy 


2ei4 


d'^Ux 

dydz 


dxdz 


+ 


d'^Uz ' 
dxdz _ 
d'^uz' 
dxdy 


-\- 2ei4 


(BIO) 

d‘^4> 

dydz' 

(Bll) 


and cyclic for Uy and Uz- Here, A is the Laplacian and e is 
the dielectric constant of the medium. For a cubic lattice, 
the mechanical boundary conditions at z = 0 explicitly 
read 


_ fduz duy\ d4>_ 

T 23 - c,^i^— + —j+ e,,— - 0, (B13) 

duz f dux duy\ 

= ■=n^+c.,(^ + ^j=0, (B14) 


while the electrical boundary condition [compare the gen¬ 
eral relation in Fq.(B9)[ leads to 


ei4 



d^\ 
dx J 



+ eokcj) 


2=0 


= 0 . 


(B15) 


In general, the wave equations can be formulated into a 
4x4 matrix A4; the condition det Ai = 0 can then used to 
find the four decay constants. In addition, the mechani¬ 
cal and electrical boundary conditions can be recast to a 
4x4 boundary condition matrix B, from which one can 
deduce the allowed phase velocities of the piezoelectric 
SAW by solving detH = 0. 

Perturbative treatment. —For materials with weak 
piezoelectric coupling (such as GaAs), the properties of 
surface acoustic waves are primarily determined by the 
elastic constants and density of the medium. Then, 
within a perturbative treatment of the piezoelectric cou¬ 
pling, one can obtain analytical expressions for the strain 
and piezoelectric fields. Here, we summarize the results 
for SAWs propagating along x'||[110] of the z||[001] sur¬ 
face following Refs. [25, 27[. Since the piezoelectric cou¬ 
pling ei 4 is small, it follows from Eq.(Bll) that cj) will be 
order 614 smaller than the mechanical displacements u, 
that is 

(l)^—u. (B16) 

e 

This results in additional terms in the wave equations 
that are of order ~ 6 ^ 4/6 « 10®N/m^. Since the elas¬ 
tic constants are 2 to 3 orders of magnitude bigger than 
this piezoelectric term, the wave equations Eqs.(BlO), 
and (cyclic versions for Uz) will be solved by the non¬ 
piezoelectric solution with corrections only at order 6 ^ 4 . 
The nonpiezoelectric solution derived in detail in Sec. A 
can be summarized as 

Ux' = 2C/Re 

Uy' — 0, 

Uz = -2iiC/Re (B17) 
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where the sound velocity v for the Rayleigh-mode follows 
from the smallest solution of 

(C44 - PV ) (CiiCii - Ci2 - CiipV ) = CiiC^p V (Cn - pV 

(B18) 

with c'li = C44 -|- (cii -I- C12) /2. The decay constant q is 
a solution of 



decay const. Q 

7 = I 7 I 

‘P 

Vs [m/s] 


GaAs 

0.50 -b 0.48i 

-0.68-b l.lOi 

1.05 

2878 

1.2 

Diamond 

0.60-b0.22i 

-1.05-b0.75i 

1.26 

11135 

0.44 


T 


Table III: Derived properties for Rayleigh surface waves, for 
both GaAs and diamond. 


(c'li - pv"^ - C44(?^) (c44 - pv"^ - ciiq‘^)+q'^ (ci2 -b 044 )^ = 0 

(B19) 

Lastly, the parameters 7 , ip can be obtained from 


7 = 


(ci2 -I- C44) q 
C44 - Ciig 2 - p-y 2 ’ 


^-2i<p 


* * 
7 - g 

1-q 


(B20) 


Now, based on the nonpiezoelectric solution given in 
Eq.(B17), the potential (p is constructed such that both 
the wave equation in Eq.(Bll) and the electrical bound¬ 
ary condition in Eq.(B15) are solved. In the {x',y',z\ 
coordinate system they read explicitly 


eA(/) 

0 


g d'^Ux' \ 

dx'dz dx'dx') ’ 


-b ei4 


dux' 

dx’ 



z=0 


(B21) 

(B22) 


One can readily check that this is achieved by the form 
proposed in Ref. [25, 27] 


2 < 0 , 


(B23) 


where (po = (ei 4 /e)C/ and Aout = i<poJ^ {0)- Here, we 
have introduced the dimensionless function (kz) which 
determines the length scale on which the electrical po¬ 
tential generated by the SAW decays into the bulk. It is 
given by 

F {kz) = 2 I All cos {/3kz + + C) + 2 ^ 36 “^^, 

(B24) 

with Ai = I All q = a + pii, and 

4 , = (B25) 

A 3 =- - —[ecostp-b eRe [Aige”*"^]-b eoRe [Aie“*‘^]] . 

e + Co 

For ARGai_ 2 ,As we obtain the following parameter val¬ 
ues [compare Ref.[25[[: |Ai| « 1.59, A 3 = —3.1, a « 
0.501, f3 « 0.472, (fi = 1.06, and ^ = —0.33. The electric 
potential for this parameter set is shown in Fig. 7. 


Appendix C: Mechanical Zero-Point Fluctuation 

In this Appendix we provide more detailed calculations 
and estimates for the mechanical zero-point motion Uq of 
a SAW. We show that they agree very well with the sim¬ 
ple estimate given in the main text. Finally, we provide 


details on the material parameters used to obtain the 
numerical estimates. 

Our first approach follows closely the one presented 
in Ref.[31[. The analysis starts out from the mechanical 
displacement operator in the Heisenberg picture 

U (x, t) [v„ (x) -b h.c.] . (Cl) 


To obtain the proper normalization of the displacement 
profiles, let us assume a single phonon Fock state, that 

is I'll) = a[j |vac) = | 0 ,..., 0 , 1 „, 0 _), where |vac) = 

njo)„ is the phonon vacuum and compute the expec¬ 
tation value of additional field energy above the vacuum 
Fniechi defined as twice the kinetic energy, since for a me¬ 
chanical mode half of the energy is kinetic, the other one 
potential [31[. We find 


Emech = 2 a;^ J (frp (r) v* (r) • v„ (r) 


= 2 / 917 max 


wr 


(C2) 

(C3) 


where the last equality defines the effective mode-volume 
for mode n. Setting Uq = max [|v„ (r)|], and assuming 
the phonon energy as Fmech = we arrive at the 

general result for a phonon mode, Ug = \/hj2pVux^', this 
confirms the simple estimate given in the main text. 

Explicit example. —Next, we provide a calculation 
based on the exact analytical results derived in Appendix 
A and B. In what follows, we assume that, in analogy to 
cavity QED, cavity confinement leads to the quantiza¬ 
tion kn = mr/Lc, where A = is the effective quan¬ 
tization area. In a full 3D model, A = LxLy where Ly 
is related to the spread of the transverse mode function 
as discussed (for example) in Refs.[14[. For simplicity, 
here we take Lx = Ly. Surface wave resonators can rou¬ 
tinely be designed to show only one resonance kg [14[. 
Within this single-mode approximation, based on results 
derived in Appendix A for a SAW traveling wave, we 
take the quantized mechanical displacement describing a 
SAW standing wave along the axis x' = (110) as 

( Xo (z) cos (kgx') \ 

0 [a -b a'l'] , (C4) 

(o(z)sin(kox') J 

Here, the functions xo {z) and /o (z) describe how the 
SAW decays into the bulk, 

Xo{z) = cos {flikgz + ip), (C5) 

Co (z) = 2 cos {Llikgz + ip + 9), (C6) 
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with material-dependent parameters = fir + 7 = 

I 7 I exp [—iff] and (/j; numerical values are presented in Ta¬ 
ble III. We note that for GaAs we find C (0) /x (0) « 1.33. 
This is in very good agreement with the numerical values 
of Ca; = \ux/4>\ = 0.98nm/V and Cz = \uz/4>\ = 1.31nm/V 
as given in Ref. [15]. Normalization of the mode-function 
allows us to determine the parameter t/g- Performing the 
integration, we find 



Uo [fm] 

so [10-®] 

00 [hV] 

Co [V/m] 

Bo [mT] 

GaAs 

1.9 

11.7 

3.1 

19.2 

— 

LiNbOa 

1.8 

11.3 

0.9-25.8 

5.8- 162.2 

— 

Quartz 

2.75 

17.3 

2.8 - 12.0 

17.3 - 75.4 

— 

Terfenol-D 

2.2 

13.8 

— 

— 

2.3 

CoFe2 04 

1.8 

11.4 

— 

— 

6.3 

Diamond 

1.17 

7.4 

— 

— 

— 


C/o = 


I h 

YTy 


(C7) 


where the parameter S depends on the material param¬ 
eters; see Table III. Using typica l mat erial parameters, 
we obtain for GaAs (diamond) ^Jflr/5 = 0.64(1.17) and 
Uq « 1.2 (1.36) fm/-y/A[/xm2]. This is in very good agree¬ 
ment with the numerical values presented in the main 
text. 

Estimates derived from literature .—In Ref.[25], it is 
shown that the SAW Rayleigh mode studied in Appendix 
A has a classical energy density £ (energy per unit sur¬ 
face area) given by 

£ = kU^H, (G 8 ) 


where U is the amplitude of the wave, k the wave vec¬ 
tor and H a material-dependent factor which is given as 
H « 28.2 X 10^°N/m^ for GaAs. By equating the clas¬ 
sical energy of the SAW given by £A, where A is the 
quantization area, with its quantum-mechanical analog 
NpiiHuj (Aph is the number of phonons), we can estimate 
the single phonon displacement Uq as 


Uo = 


hu! 

kHA 


= 1,05X10-^, (C9) 


with U = Uq sJ Aph. This estimate is also found to be in 
very good agreement with a result given in Ref. [38] as 


Un = C, 


I 2h 

pvsA 


1.7x10-21 (GIO) 

VA 


where C is a normalization constant with numerical value 
C « 0.45 for GaAs [38] . Therefore, for an effective mode 
area of = '/A = l/rm we find a single phonon displace¬ 
ment of Uq « Ifm. This confirms the estimates given in 
the main text. 


Table IV: Estimates for zero-point fluctuations (mechanical 
amplitude Uq, strain sq , electrical potential (f)Q, electric field 
^0 and magentic field Bq) close to the surface {d <C A) for typ¬ 
ical piezo-electric and piezo-magnetic (magnetostrictive) ma¬ 
terials. All values must be multiplied by the universal scaling 
factor XjA [/rm^j; thus, they refer to an effective surface 
mode area of size A = Lower (upper) bounds for </)o 

and ^0 comprise minimum (maximum) non-zero element of e 
with maximum (minimum) non-zero element of e. We have 
set k = 2tx j [xm. Details on cut-direct ions and material pa¬ 
rameters are given in the text. 


correspondence to Eqs.(2) and (3), with the appropriate 
replacements (using standard notation) E —^ H, D —> B, 
^ij —t Pij and Cijk —>■ hijk [32]. Goupling between me¬ 
chanical and magnetic degrees of freedom is described 
by the piezomagnetic tensor h which can reach values as 
high as ^ 700T/strain [35] ; for our estimates we have 
referred to Terfenol-D, where /115 « 167T/strain. The 
magnetic field associated with a single phonon can then 
be estimated as Bq « higSg, where hiQ refers to a typical 
(non-zero) element of h. 

For the piezoelectric materials GaAs, LiNbOg and 
Quartz all material parameters have been obtained from 
Ref. [26] . Phase velocities for typical cut directions have 
been used, that is (100) [001] GaAs, Y-Z LiNbOg and 
ST Quartz. For the piezomagnetic (magnetostrictive) 
materials CoFe 204 and Terfenol-D all material param¬ 
eters have been taken from Ref. [34] . We have used the 
phase velocities of the bulk shear waves given in there 
as Ush = 3.02 X lO^m/s and Ush = 1-19 x lO^m/s for 
CoFe 2 04 and Terfenol-D, respectively. This gives an 
conservative estimate for Uq, since Rayleigh modes have 
phase velocities that are lower than the ones of bulk 
modes [14] . For example, in the case of CoFe 2 04 in 
Ref. [35] wave velocities for Rayleigh-type surface waves in 
a piezoelectric-piezomagnetic layered half space are found 
to be Vs « 2840m/s < Ush- 


Appendix D: Zero-Point Estimates 

In this Appendix we provide details on piezo-magnetic 
materials and numerical estimates of the zero-point quan¬ 
tities for several relevant materials. The results are sum¬ 
marized in Tab. IV. The underlying input parameters are 
given below. 

Theoretically, piezo-magnetic materials with a large 
magneto-strictive effect are typically described in a 1:1 


Appendix E: SAW Cavities 

In this Appendix, we present a detailed discussion of 
the theoretical model describing the SAW resonator. 

Typically, a SAW cavity is based on an on-chip dis¬ 
tributed Bragg reflector formed by a periodic array of 
either metal electrodes or grooves etched into the sur¬ 
face; see Fig.l. In such a grating, each strip reflects only 
weakly, but, for many strips A ^ 1, the total reffec- 
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groove depth hjX c [10 ^ 

Figure 8: (color online). Total reflection coefflcient |i?| as a 
function of the normalized groove-depth /i/Ac for N = 100 
(blue dashed) and N = 300 (red solid). Here, vjjp = 0.5 and 
material parameters for LiNbOa have been used (see text). 

tion |i?| can approach unity if the pitch p equals half the 
wavelength, p = Ac/2. This Bragg condition defines the 
center frequency 

fc = Vs/2p. (El) 

At / = fc, the total reflection coefficient is given by 

|i?| = tanh(iV |rs|), (E2) 

where N is the number of strips and is the reflec¬ 
tion coefficient associated with a single strip [14, 15]. 
The total reflection coefficient |i?| goes to unity in the 
limit iVjrsI 1; see Fig. (8). Typically, N > 200 
and [r^l « (1-2)% [14]. For / « fc, jr^j in¬ 
creases with the normalized groove depth as [rg] = 
Cih/XcSin {Trw/p)+C2 (h/\c) cos (irw/p), with material- 
dependent pre-factors [37] . For LiNbOa, Ci = 0.67 and 
C 2 = 42 [37]. As argued in Ref. [37], the first term ^ Ci 
is due to a impedance mismatch, while the second one 
~ C 2 is due to the stored energy effect. 

Due to the distributed nature of the mirror, strong re¬ 
flection occurs over a fractional bandwidth only, given 
by bf/fc ~ 2 |rs|/ 7 r. In practice, the cavity formed 
by two reflective gratings can be viewed as an acous¬ 
tic Fabry-Perot resonator with effective reflection cen¬ 
ters, sketched by localized mirrors in Fig.l, situated 
at some effective penetration distance into the grating, 
given by Lp = tanh [(A^ - 1) [rg]] Ac/(4 jr^j) « Ac/4|rgj 
[14, 15, 36] . Therefore, the total effective cavity size 
along the mirror axis is Lc ^ D + 2Lp, where D is the 
physical gap between the gratings; compare Fig.l. For 
N « 100—300, h/Xc « 2%, we then obtain Lc « 38Ac and 
Lc « 42Ac for D — 0.75Ac and D = 5.25Ac, respectively. 
In analogy to an optical Fabry-Perot resonator, the mode 
spacing can then by estimated as Af /fc = Xc/2Lc ~ jrgj. 


Since this is larger than 6f/fc, SAW resonators can be 
designed to host a single resonance only [14] . 

The total decay rate of this resonance k can be decom¬ 
posed into four relevant contributions [37], k = Kbfe + K£; + 
+ which inludes conversion into bulk modes ^ Kbk, 
diffraction losses ^ Kd, internal losses due to material im¬ 
perfections ^ Km, and leakage (radiation) losses due to 
imperfect mirrors ~ The associated Q-factors given 
by Qi = tJc/ni- The desired decay rate is Kgd = Kr, 
whereas the undesired one is Kbd = 

Here, Kd is associated with diffraction losses due to spill¬ 
over beyond the aperture of the reflector. It can be 
made negligible by lateral confinement using for exam¬ 
ple waveguide structures, focusing or etching techniques 
[14, 38, 39] . Qm refers to losses due to interaction with 
thermal phonons, losses due to defects in the material 
and propagation losses due to contamination [36, 37] . 
These losses ultimately limit Q: Low temperature experi¬ 
ments on quartz have demonstrated SAW resonators with 
QmXf [GHz] > 10® [21, 22] . Another source of losses is 
due to mode-conversion into bulk-modes. Measurements 
show that Qbk = 27riVeff/[C'b (/i/Ac)^] with Wff = Lc/Xc 
and a material-dependent pre-factor Ct [85] ; for LiNbOa 
(Quartz), Cb = 8.7(10), respectively [85, 86 [. Typically, 
Kbk is found to be negligible for small groove depths, 
h/Xc < 2% [37]. Finally, Kr arises from leakage through 
imperfectly reflecting gratings (|i?| < 1 ); in direct anal¬ 
ogy to optical Fabry-Perot resonators, the associated Q- 
factor is given by Qr = 27rA^eff/ (l — l-Rp)- Assuming 
negligible diffraction losses (that can be minimized via 
waveguide-like confinement [14, 42[) and cryostat tem¬ 
peratures, the total Q-factor is then given by 

Q ^=Qm+Qf,i}+Qr^- (E3) 


Appendix F: Charge Qubit Coupled to SAW Cavity 
Mode 

We consider a GaAs charge qubit embedded in a 
tunnel-coupled double quantum dot (DQD) containing 
a single electron. In the one-electron regime the single¬ 
particle orbital level spacing is on the order of ~ ImeV. 
Therefore, the system is well described by an effective 
two-level system: The state of the qubit is set by the po¬ 
sition of the electron in the double-well potential, with 
the logical basis \L) , |i?) corresponding to the electron 
localized in the left (right) orbital. The Hamiltonian de¬ 
scribing this system reads 

Hch = + tca^, (FI) 

with the (orbital) Pauli-operators defined as ct* = 
|L) {L\ — |i?) (i?| and ct^ = |L) {R\ + |i?) {L\, respectively. 
In Eq. (5) , e refers to the level detuning between the dots, 
while tc gives the tunnel coupling. The level spli tting be- 
tween the eigenstates of TJch is given by H 
with a pure tunnel-splitting of Q = 2tc at the charge 
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degeneracy point (e = 0 ); typical parameter values are 
tc ^ /reV and e ~ /reV, such that the level splitting 
n ~ GHz lies in the microwave regime. At the charge 
degeneracy point, where to first order the qubit is insen¬ 
sitive to charge fluctuations {dfl/de = 0 ), the coherence 
time has been found to be T 2 « 10ns [45]. 

We now consider a charge qubit as described above in¬ 
side a SAW resonator with a single resonance frequency 
ojc close to the qubit’s transition frequency, uJc ~ fl, that 
is the regime of small detuning 5 = — Wc ~ 0; note 

that single-resonance SAW cavities can be realized rou¬ 
tinely with today’s standard techniques [14]. Within this 
single-mode approximation, the Hamiltonian describing 
the SAW cavity simply reads 

Tfcav — (F 2 ) 

where (a) creates (annihilates) a phonon inside the 
cavity. The electrostatic potential associated with this 
mode is given by (p (x) = p (x) [a -I- a^] , where the mode- 
function (j) (x) can be obtained from the correspond¬ 
ing mechanical mode-function w (x) via the relation 
eA(/)(x) = CkijdjdkWi (x); here, A is the Laplacian, Ckij 
the piezoelectric tensor and e the permittivity of the ma¬ 
terial. The electron’s charge e couples to the phonon in¬ 
duced electrical potential </>. In second quantization, the 
piezoelectric interaction reads Hint — dx0(x)h(x), 

where e is the electron’s charge, n (x) = ipl. (x) tpa (x) 
is the electron number density operator and (x) cre¬ 
ates an electron with spin a at position x [24] . Since p (x) 
varies on a micron length-scale which is large compared to 
the spatial extension ^ 40nm of the electron’s wavefunc- 
tion in a QD [64], the electron density is approximately 
given by a delta-function at the center of the correspond¬ 
ing dots. For the DQD system under consideration Hint is 
then approximately given by Hint = e 4> (^i) here, 
Xi refers to the center of the electronic orbital wave- 
function ipi (x) of dot i = L,R. Note that this form 
of Hint becomes exact if the overlap integral vanishes, 
that is if / dx(j) (x) tjj’l (x)'ipfi (x) = 0 is satisfied. As 
shown below, for a mode-function p (x) of sine-form, this 
condition maximizes the piezoelectric coupling strength 
between the electronic DQD system and the phonon 
mode. For the charge qubit system under considera¬ 
tion coupling to the cavity mode is then described by 
Hint = e[a +aP)[(j){xL)\L) {L\+(j){xR)\R) {R\]] here, 
Xi refers to the center of the electronic orbital wavefunc- 
tion ipi (x) of dot i = L,R and the transverse direction 
y has been integrated out already. To obtain strong cou¬ 
pling between the qubit and the cavity, we assume a mode 
profile 4>{x) = (/?osin (fca;), with a node tuned between 
the two dots, such that p (xl) = <posin {kl/2) = —p (xr); 
here, I gives the distance between the two dots. Note that 
the single phonon amplitude, defined as ipo = cptjR (kd), 
with R {kd) « 0 for d ^ A, accounts for the decay 
of the SAW into the bulk. For A = (0.5 — 1) /im and 
a 2DEG (where the DQD is embedded) situated a dis¬ 
tance d = 50nm below the surface, however, the single 


phonon amplitude is reduced by a factor of ^ 2 only, 
ipo « (0.45 — 0.52) (po; see Appendix B for details. Then, 
the coupling between qubit and cavity reads 

Hint = Och {a + a'^) ^ , (F3) 

where the single-phonon coupling strength is 

gch = epoR {kd) sin {kl/2) « ^===- (F4) 

[/rm2] 

Here, we have assumed A « 2/ such that the geomet¬ 
rical factor sin( 7 r//A) « 1 [46]. In principle, the cou¬ 
pling strength ^ch could be further enhanced by addi¬ 
tionally depositing a strongly piezoelectric material such 
as LiNbOa on the GaAs substrate. Moreover, compari¬ 
son with standard literature shows that the piezoelectric 
electron-phonon coupling strength can be expressed as 
Spe = VPUq « e(ei 4 /e)f 7 o = epo, where P = {eeu/e)'^ 
is a material parameter quantifying the piezoelectric cou¬ 
pling strength in zinc-blend structures [87, 88 ] . Using 
P = 5.4 X for GaAs, the single phonon Rabi 

frequency can be estimated as gpe « 2.87geV/^/A[yuP]. 
This corroborates our estimate for gnt- 

In summary, the total system can be described by the 
Hamiltonian H = Hch + 77cav + Hint, 

H = -I- tccr^ + WcO^a -|- gch (a + G cr^. (F5) 

This corresponds to the generic Hamiltonian for a qubit- 
resonator system [31] . It is instructive to rewrite H in 
the eigenbasis of Hch, given by 

|-l-) = siii9\L)+cos9\R), (F 6 ) 

|—) = cos9\L) — sm9\R), (F7) 

where the mixing angle 9 is defined via tan 9 = 
2A/(e + U). In a rotating wave approximation 
{h,gch ‘C Wc), H then reduces to the well-known Hamil¬ 
tonian of Jaynes-Gummings form 

Of 

H^6S^+gch^{S+a + S-a^), (F 8 ) 

where 5 = D — Wc, = (|-|-) (-|-| — |—) (—|) /2 and = 

l±) (T|. 

Appendix G: SAW-Based Cavity QED with Spin 
Qubits in Double Quantum Dots 

In this Appendix, we show in detail how to realize the 
prototypical Jaynes-Gummings dynamics based on a spin 
qubit encoded in a double quantum dot (DQD) inside a 
SAW resonator. 

We consider a double quantum dot (in the two-electron 
regime) coupled to the electrostatic potential generated 
by a SAW. The system is described by the Hamiltonian 

HdQD — Hq Hcnv T Hint, 


(Gl) 
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where -ffcav and i/i„t describe the DQD, the cav¬ 
ity and the electrostatically mediated conpling between 
them, respectively. In the following, the different contri- 
bntions are discnssed in detail. 

Double Quantum Dot .—The DQD is modeled by the 
standard Hamiltonian 

Ho = Hc + Ht+Hz. (G2) 

Here, He gives the electrostatic energy 


\T+) Sf,. = +1 



He = '^ e-i^ia + U + ULRULnn, (G3) 

1,(7 i = L,R 

where (due to strong confinement) both the left and right 
dot are assumed to support a single orbital level with en¬ 
ergy Ci (z = L, R) only; U and Ulr refer to the on-site 
and interdot Goulomb repulsion, respectively. As usual, 
nia = d\^dia and Ui = +nii refer to the spin-resolved 
and total electron number operators, respectively, with 
the fermionic creation (annihilation) operators c?|^ (dia) 
creating (annihilating) an electron with spin cr =t) i iii 
the orbital i = L,R. We focus on a setting where an 
applied bias between the two dots approximately com¬ 
pensates the Goulomb energy of two electrons occupy¬ 
ing the right dot; that is, €l « en + U — Urr- Thus, 
we restrict ourselves to a region in the stability diagram 
where only (1,1) and (0,2) charge states are of interest. 
All levels with (1,1) charge configuration have an elec¬ 
trostatic energy of + ^r + Urr, while the 

(0, 2) configuration has A(o, 2 ) = 2e_R -|- U. As usual, we 
introduce the detuning parameter e = tL — ^R + URR — U. 
In this regime, the relevant electronic levels are defined 
as |T+) = IM), |T_) = IM), |To) = (IM) + IM))/V2, 
I'S'ii) = (IM) - |M))/v^ and IS' 02 ) = d^j^.^d^j^\0) with 

io)- 

Next, Ht describes coherent, spin-preserving interdot 
tunneling 

Ht = tc^d\^dRa +h..c., (G4) 

<7 

where tc is the interdot tunneling amplitude. Lastly, Hz 
accounts for the Zeeman energies, 

Hz = 9t^B Y. 

i—L,R 

where g is the electron g-factor and gtR the Bohr 
magneton, respectively. In the presence of a micro- 
/nanomagnet, the two local magnetic fields are in¬ 
homogeneous, Bi 7 ^ Br. We can then write B^ = 
Bq + Bm (xi), where Bg is the external homogeneous 
magnetic field, while B^ (x^) is the micromagnet slant¬ 
ing field at the location of dot x^. In practice, Bq is a 
few Tesla, at least larger than the saturation field of the 
micromagnet Bq > 0.5T, while the magnetic gradient 
AB = |jB„i (xfl;) — Bm (x 7 ^)|| can reach AB « lOOmT, 
corresponding to an electronic energy scale of \gfj.BAB\ « 


Figure 9: (color online). Illustration of the relevant electronic 
levels under consideration. The triplet levels with Stot 7^ 1 
can be tuned off-resonance by applying a sufficiently large 
homogeneous magnetic field. 

2/ieV [89]. Field derivatives realized experimentally are 
dBjn.zIdx « 1.5mT/nm. Alternatively, the magnetic 
gradient can be realized via the Overhauser field, as ex¬ 
perimentally demonstrated for example in Ref. [ 66 ] . 

Note that the Fermi contact hyperfine interaction be¬ 
tween electron and nuclear spins reads H^f = hi • S^. 
Here, hi is the Overhauser field in QD i = L,R. When 
treating hi as a classical (random) variable, Huf is equiv¬ 
alent to Hz and thus one can absorb hi into the definition 
of the magnetic field Bi in Eq.(G5); also see Ref. [89]. 

To facilitate the discussion, we introduce the mag¬ 
netic sum field B = (B^ -I-B^j) /2 and the difference 
field AB = (Bfl - Bi) /2. While B conserves the total 
spin, that is [B (Si-b S^^), (SlS^)^] = 0, the gra¬ 
dient field AB does not. We set the quantization axis z 
along B = Bz. For sufficiently large magnetic field B the 
electronic levels with = 5'£ -b yf 0 are far detuned 
and can be neglected for the remainder of the discussion. 
Therefore, in the following, we restrict ourselves to the 
= 0 subspace. The components AB“’^ give rise to 
transitions out of the (logical) subspace = 0. Since 
these processes are assumed to be far off-resonance, they 
are neglected leaving us with the only relevant magnetic 
gradient A = g^BAB^/2; compare also Refs. [65, 66 [. 
For a schematic illustration, compare Fig. 9. 

In summary, in the regime of interest Hq simplifies to 

Ho = |(|^o2)(5ii|-bh.c.)-e|5o2)(5o2| (G 6 ) 

-A(|To)(^n|+h.c.). 

The eigenstates of Hq within the relevant = S'£ -b 
S'£. = 0 subspace can be expressed as 

\l) = ai\To)+Pi\Sn) + Ki\So2), (G7) 

with corresponding eigenenergies e/(Z = 0,1,2). The 
spectrum is displayed in Fig. 10. For large negative de¬ 
tuning —e ^ tc, the level | 2 ) is far detuned, and the 
electronic subsystem can be simplified to an effective two- 
level system comprising the levels {| 0 ) , | 1 )}, that is 

ffo«y(|l)(l|-|0)(0|), (G8) 
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Figure 10: (color online). Spectrum of Hq for ta = 5A. The 
three eigenstates |Z) are displayed in green dotted {1 — 2), 
red solid (Z = 1) and blue dashed (Z = 0), respectively. The 
triplets |T±) are assumed to be far detuned by a large exter¬ 
nal field and not shown. For large negative detuning e > —tc, 
the hybridized levels {|0) , |1)} can be used as qubit. 

which can be identified with a ’singlet-triplet’-like logical 
qubit subspace. Here, loq = ei — eg refers to the qubit’s 
transition frequency. Note that the magnetic gradient 
causes efficient mixing between ITq) and IIS'!!) for A > 
IZg/el. In the regime of interest, the dominant character 
of the qubit’s levels is |1) « lll'fl'), |0) « I'fl'.y.) (or vice 
versa) [65] and the transition frequency is approximately 
ojo ~ 2A. For A « l^eV, the transition frequency ojq = 
ei — eo ~ 2/reV « 3GHz matches typical SAW frequencies 
~ GHz. 

Coupling to SAW phonon mode .—Along the lines of 
Appendix F, again we consider a SAW resonator with 
a single relevant confined phonon mode of frequency Wc 
close to the qubit’s transition frequency uig. For a DQD in 
the Zwo-electron regime, in the basis of Eq. (G6) coupling 
to the resoantor mode can be written as [46] 

Hint = 5o [a + G l‘S' 02 ) (>S'o 2 | j (G9) 

where gg = eipgijggn- Here, rjggo is a geometrical fac¬ 
tor accounting for the DQD’s position with respect 
to the mode-function (/>(x); it is defined according to 
‘PoVgeo = (f> (^r) ~ For example, taking a stand¬ 

ing wave pattern along x as demonstrated experimentally 
in Ref. [90], together with a transverse mode function re¬ 
stricting the spread in the y-direction [39, 42], we obtain 
rjgeo = sin (27ra:i?,/A) — sin (27ra;L/A). It takes on its maxi¬ 
mum value ?7opt, when tuning a node of the standing wave 
at the center between the two dots, that is xr = 1/2, 
XL = —^/2; this gives ryopt = 2sin(7rZ/A), where I is the 
distance between the two dots [46] ; as compared to the 
charge qubit described in Appendix F, there is an addi¬ 
tional factor of two, since here we consider a DQD in the 
Zwo-electron regime, whereas the charge qubit consists of 
one electron only. For typical parameters (Z = 220nm, 
A « 1.4/rm) as used in Ref.[46[, we get rjopt ~ 0.95, while 
I = 220nm, A « 0.5yxm leads to the largest possible value 


Figure 11: (color online). The product kqki directly affects 
the effective single-phonon Rabi frequency gQo/go ~ kqKi 
[ 69] , while the difference | fti — kq | determines the robustness 
of the qubit against charge noise. Here, tc = 5A. 

Vopt ~ 2. 

In summary, within the effective electronic two-level 
subspace {|0), |1)}, the system is described by the Hamil¬ 
tonian 

Huqd = E ("' + 1^) (^1 (GIO) 

+ KgKiVpe (|0) (1| -|- h.C.) -|- UJcO^a, 

where Vpe = 50 [a + ■ Applying a unitary trans¬ 

formation to a frame rotating at the cavity frequency 
coc according to Hdqd = UHbqdU^ + iUU\ with 
U = exp [iuJct (a'^a -I- ^ | 1 ) ( 1 | — ^ |0) (0|)], performing 
a rotating-wave approximation (RWA), and dropping 
a global energy shift e = (eo + ei)/2, we arrive at 
the effective (time-independent) Hamiltonian of Jaynes- 
Cummings form 

Hbqd = SS^+gQu[S+a + S-a^, (Gil) 

where we have introduced the spin operators S~^ = | 1 ) (0| 
and = (|1) (1| — |0) (0|) /2. Moreover, S = ujg — 
is the detuning between the qubit’s transition frequency 
ojg and the cavity frequency Wc, and the effective single¬ 
phonon Rabi frequency is defined as 

5qd = KoR-iVgeoec/giF {kd) « 2KoKi5ch- (G12) 

The coupling between the qubit and the cavity mode is 
mediated by the piezoelectric potential; therefore, it is 
proportional to the electron’s charge e and the single¬ 
phonon electric potential (j)g. Due to the prolonged deco¬ 
herence timescales, here we consider an effective (singlet- 
triplet like) spin-qubit rather than a charge qubit, such 
that the coupling ^qd is reduced by the (small) admix¬ 
tures with the localized singlet ki = (5'o2|Z). Increas¬ 
ing KgKi leads to a stronger Rabi frequency gQu, but 
an increased difference in charge configuration | — Kq | 

makes the qubit more susceptible to charge noise. For 
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cavity size ^[pm^] 

PQD [MHz] 

Q 

c 

small cavity 

1 

200 

10 

4.25 

large cavity 

500 

9 

lO'^ 

8.5 


Table V: Estimates of the single spin cooperativity C for a 
DQD singlet-triplet qubit with « 100ns, in a SAW cavity 
at gigahertz frequencies a;c/27r « 1.5GHz and cryostat tem¬ 
peratures where hth ~ 0 for both a small, low Q and large, 
high Q SAW-resonator. The coupling strength qqu could be 
further increased by additionally depositing a strongly piezo¬ 
electric material such as LiNbOs on the GaAs substrate and 
spin-echo (and/or narrowing) techniques allow for dephasing 
times extended by up to three orders of magnitude [63, 68]. 


typical numbers ~ 5/ieV, e « — 7/ieV, A « l/reV), we 
get kqKi « 5 X 10“^, \k\ — KqI ~ 2 x 10“^; see Fig. 11. 
For I « 250nm, A « 0.5/im, and d « 50nm we can then 
estimate 


ffQD/S 


200MHz 

\/H[/rm2] 


(G13) 


We take this coupling strength as a conservative estimate, 
since optimization against the relevant noise sources as 
done in Ref.[69] yields an optimal point with kqKi « 
0.3 and « 1.5GHz. Resonance (h = 0) yields a 

SAW wavelength of A « 2^m; accordingly, for a fixed 
dot-to-dot distance I = 250nm, « 0.76 (whereas 

a larger DQD size I = 400nm as used in Ref. [24] gives 
Vgeo ~ 1.18). For I = 250nm, d = 50nm, we then obtain 
gQB/h « GOOMHz/yif^ m2], which is a factor of three 
larger than the estimate quoted above. 

Cooperativity .—In this context, an important figure 
of merit is the single spin cooperativity [60], C = 
5 qd^ 2 /k (h-th + 1), where k, = uJc/Q is the mechanical 
damping rate and hth = \ -g equilib¬ 

rium phonon occupation number at temperature T; here, 
since hwc 3> ksT for cryostatic temperatures, hth ~ 0. 
For singlet-triplet qubits in lateral QDs, ~ 100ns [63]; 
using spin-echo techniques, experimentally this has even 
been extended to T 2 = 276/rs. Even in the absence of 
spin-echo pulses, with a far-from-optimistic dephasing 
time T 2 « 100ns [68], for a moderately small cavity size 
A « 100/im^, a quality factor of Q = 900 is sufficient to 
reach C « gQj^T^QlLVc « 3.8. Note that C > 1 allows to 
perform a quantum gate between two spins mediated by 
a thermal mechanical mode [9]. 

Discussion of approximations .—To arrive at the effec¬ 
tive Hamiltonian given in Eq.(Gll), we have made two 
essential approximations: (i) first, we have neglected the 
electronic level |2) yielding an effective two-level system 
(TLS), and (ii) second, we have applied a RWA leading 
to a major simplification of FfDQDi see Eq.(Gll) as com¬ 
pared to Eq.(GlO). In order to corroborate these approx¬ 
imations, we now compare exact numerical simulations of 
the full system where none of the approximations have 
been applied to the simplified, approximate description 
described above. While the dynamics of the former is 



Figure 12: (color online). Exact numerical simulations of the 
full (blue, solid line) and approximate (cyan circles) Master 
equations as given in Eqs.(G14) and (G15) respectively. Plots 
are shown for (a) the electronic inversion (5^)^, (b) the cavity 
occupation (n)j and (c) the error Tr [p |2) (2|] quantifying the 
leakage to |2). The latter is found to be negligibly small ~ 
10“®. We have set <5 = wo — Wc = 0. Numerical parameters; 
tc = 10/reV, e = — 7/ieV, A = 1/reV, Pgeoetpo = 5.2 x 10“^/ieV 
such that gQD = 4 X 10“^^eV ~ 6 MHz. The cavity decay rate 
is K = 5 rQD/ 2 , corresponding to Q ~ 10^. 


described by the Master equation 

p = -i [Hq + FIcav + 77int, p] + CD [a] p, (G14) 

with TJo, FIcav and Flint given in Eqs.(G6), (F2) and (G9), 
respectively, the latter is described by a similar Master 
equation with the coherent Hamiltonian term replaced 
by the prototypical Jaynes-Gummings Hamiltonian, 

p= -i [55'"' -I- pQD (5+a -I- S~a^) , p] + kV [a] p, (G15) 

with p referring to the density matrix of the combined 
system comprising the DQD and the cavity mode. Here, 
we have also accounted for decay of cavity phonons out 
of the resonator with a rate k, described by the Lindblad 
term D[a\p = apa) — ^ {a^a, p}. As a figure of merit 
to validate approximation (i) we determine the popula¬ 
tion of the level level |2), that is Tr [p |2) (2|], describing 
the undesired leakage out of the logical subspace; ide¬ 
ally, this should be zero. Note that leakage into the 
triplet levels \T±) could be accounted for along the lines, 
but they can be tuned far off-resonance by another, in¬ 
dependent experimental knob, the external homogenous 
magnetic field. The results are summarized in Fig. 12: 
We find very good agreement between the exact and 
the approximate model, with a negligibly small error 
Tr[p|2)(2|] ~ 0(10“®). This justifies the approxima¬ 
tions made above and shows that (in the regime of inter¬ 
est) the system can simply be described by Eq.(G15). 
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Figure 13: (color online). Transition frequency cjq (blue solid) 
and its sensitivity against charge noise induced fluctuations 
in e for intermediate (a) and large (b) negative detuning; for 
large negative detuning luq ~ 2A, and the sensitivity duJo/de 
practically vanishes leaving nuclear noise as the dominant de¬ 
phasing process. By occasional refocusing of the spin states, 
this regime can be used for long-term storage of quantum 
information [69]. Other numerical parameters: tc = 5A. 


in the Overhauser field generated by the nuclear spins 
[65, 66, 68]. In our model, this can directly be iden¬ 
tified with a random, slowly time-dependent parameter 
(5A = SA {t). In the relevant subspace {|0), |1)}, SH^uc 
is given by 6 Hnuc = -SAY,k,i <^kl3i [\k) {l\ + h.c.]. Typi¬ 
cally, SA « 0.1/xeV ujc is fulfilled, such that we can ap¬ 
ply a RWA yielding SH^nc « -‘^SAY.i ail^i 10 (0; Phys¬ 
ically, SHnuc is too weak to drive transitions between 
the electronic levels |0) and |1) which are energetically 
separated by wn ~ uJr- Then, in the spin basis used in 
Eq.(Gll), we find 


6 Hnnc=Sit)S\ (G18) 

where the gradient noise is given by S (t) = 

2SA{t) {aoPo - OiiPi). For (tc,e. A) « (5,-7,1)/xeV, 
ao/3o — ctiPi ~ 0.9. Therefore, when also accounting for 
nuclear noise as described by Eq.(G18), the full model 
[compare Eq.(Gll)] reads 

Hbqb = [5 + <5 (t)] 5" + ffQD [S+a + 5-ai] . (G19) 


Noise Sources for the DQD-based System 

Charge noise .—In a DQD device background charge 
fluctuations and noise in the gate voltages may cause un¬ 
desired dephasing processes. In a recent experimental 
study [91], voltage fluctuations in the intedot detuning 
parameter e have been identified as the main source of 
charge noise in a singlet-triplet qubit. Gharge noise can 
be treated by introducing a Gaussian distribution in e, 
with a variance typically CTe « (1 — 3) fieV [89]. The 
qubit’s transition frequency wq, however, turns out to 
be rather insensitive to fluctuations in e, with a (tun¬ 
able) sensitivity of approximately dcoo/de < 10“^; see 
Fig. 13. In agreement with experimental results presented 
in Ref. [91], we find dojo/de ~ wq, indicating cjq to be an 
exponential function of e. At very negative detuning e, 
dephasing due to charge noise is practically absent, and 
T 2 will be limited by nuclear noise [91]. Fluctuations 
in the tunneling amplitude tc can be treated along the 
lines: we find wn to be similarly insensitive to noise in 
dujo/dtc - 10 - 2 . 

Nuclear noise: Spin echo .—The electronic qubit in¬ 
troduced above has been defined for a fixed set of pa¬ 
rameters (tc,e, A); compare Eq.(G6). Now, let us con¬ 
sider the effect of deviations from this fixed parame¬ 
ters, Hq — Hq + SH, where SH can be decomposed as 
SH = SHei + SHciuc with 

SHci = ^ ( 1 ^ 02 ) (5ii|+h.c.)-5e|5o2)(^02f(?16) 

SHnuc = -<5A(|To)(^ii|+h.c.), (G17) 

where Stc and Se can be tuned electrostatically and ba¬ 
sically in-situ. In most practical situations this does not 
hold for i5A: The primary source of decoherence in this 
system has been found to come from (slow) fluctuations 


Since the nuclear spins evolve on timescales much longer 
than all other relevant timescales ^ [64], the 

Overhauser noise term can be approximated as quasi¬ 
static, that is S (t) = S. As experimentally demonstrated 
in Ref. [66], the slow (nuclear) noise term ~ <5 (t) can be 
neutralized by Hahn-echo techniques. Here, the dephas¬ 
ing time of the electron spin qubit was extended by more 
than three orders of magnitude from ~ (10 — 100) ns 
to T 2 = 276/rs. 

In the following, assuming nominal resonance <5 = 0, 
we detail a sequence of Hahn-echo pulses that cancels 
the undesired noise term and restores the pure, res¬ 
onant Jaynes-Gummings dynamics: We consider four 
short time intervals of length r, for which the unitary 
evolution is approximately given by [7^ « 1 — iHiT, in¬ 
terspersed by three 7r-pulses. First, we let the system 
evolve with Hi = SS^ + gqu [S''''a -|- S'-af], then we ap¬ 
ply a TT-pulse along x (S^ S^, Sy -Sy) 

such that —>■ and the system evolves 

in the second time interval with H 2 = —dS^ + 
9qd [<S'-a -I- 5'+ of]. Next, we apply a 7r-pulse along z 
{S^ ^ SC ^ -S‘^,sy -sy) such that S'± ^ -S'± 
leading to = —SS^ — (/qd [5-a-I-S'+a^]. Finally, 
a TT-pulse along y {S^ -S^S^ -S^,Sy Sy) is 

applied such that S^ —>■ —S^ giving H 4 = -I- 

9 QD [5'“'"a -I- S'-ai]. In summary, the system evolves over 
a time interval of 4r according to t/eff = U 4 IJ 3 U 2 U 1 « 
1 — ir Hi = ^ — iHeffir with the effective Hamilto¬ 
nian 

i?eff = ^ [5+a + 5-ai] . (G20) 

Thus, in order to cancel the noise term, the effective 
single-phonon coupling strength is only lowered by a fac¬ 
tor of 1/2. 
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Different Spin-Resonator Coupling 

In Appendix G, we have shown how to realize the 
prototypical Jaynes-Cummings Hamiltonian for SAW 
phonons interacting with a DQD; see Eq.(Gll). Alterna¬ 
tively, if one does not absorb the gradient A into the defi¬ 
nition of the qubit basis, one can identify the logical sub¬ 
space with the electronic states |To) and [S'), where \S) 
is one of the two hybridized singlets (while the other one 
|S") is far detuned and neglected) [65, 66, 68, 91]. Here, 
the electronic Hamiltonian reads Hq = —J{e)\S){S\ — 
A (|To) (S'!-I-h.c.), where A = (S'ii|S')A and J (e) de¬ 
scribes the exchange interaction. In this regime, the spin- 
resonator interaction takes on a form that is well known 
from other (localized) implementations of mechanical res¬ 
onators [9], namely 

Hint = gqd {a + a^) \S) {S\, (G21) 

which can be viewed as a phonon-state dependent force, 
leading to a shift of the qubit’s transition frequency de¬ 
pending on the position x = (a-I-a^) /v^- Here, the 
single-phonon Rabi frequency is = K^r]gene(j)Qj- (kd), 
with Ks = {So 2 \S). Based on the coupling of Eq.(G21), 
one can envisage a variety of experiments known from 
quantum optics: For example, in the limit of vanishing 
gradient A = 0, the x quadrature of the phonon mode 
could serve as a quantum nondemolition variable, as it 
is a integral of motion of the coupled system of phonon 
mode and electronic meter. 


Appendix H: Generalized Definition of the 
Cooperativity Parameter 

In this Appendix we provide a generalized discussion 
of the cooperativity parameter C which in particular ac¬ 
counts for losses of the cavity mode other than leakage 
through the non-perfect mirrors. Furthermore, we de¬ 
rive a simple, analytical estimate for the state transfer 
fidelity JF in terms of the parameter C and undesired 
phonon losses with a rate ~ Kbd- 

We consider a single qubit {10), 11)} coupled to a cavity 
mode. The system is described by the Master equation 

g = {Kgd + Kbd)T> [a] g-i [iJjc, £»] + rdephT> [|1) (1|] g, 

Cog Vg 

(HI) 

where 'D[a\g = aga^ — ^ {o^a, £>}. The first term de¬ 
scribes decay of the cavity mode. The corresponding de¬ 
cay rate can be decomposed into desired (leakage through 
the mirrors) and undesired (bulk mode conversion etc.) 
contributions, labeled as Kgd and Kbdi respectively. Thus, 
we write k = Kgd -f Kbd- The second term with (on reso¬ 
nance) Hdc = g (5'+a -|- S'“a^) refers to the coherent in¬ 
teraction between qubit and cavity mode, while the last 
term describes pure dephasing of the qubit with a rate 

Edeph- 


In the bad cavity limit (where k 3> g,rdeph)> one can 
adiabatically eliminate the cavity mode by projecting the 
system onto the cavity vacuum, Vg = Tr^av [f?] G Pcav = 
p® |vac) (vac|. Standard techniques (perturbation theory 
up to second order in V, compare Ref. [92]) then yield the 
effective Master equation for the qubit’s density matrix 
P = Trcav [Vg] only, 

p = kV[S-]p + T^,pt,V[\l){l\]p, (H2) 

with the effective decay rate k = Ag'^ Jk. 

For comparison, the same procedure in standard cavity 
QED, where k = Kgd and TdephT^ [|1) (1|] p jV [S'”] p, 
yields the effective Master equation for the atom only, 
p = kV [S”] p + jV [S”] p. Therefore, the atom decays 
with an effective spontaneous emission rate ytot enhanced 
by the Purcell factor, ytot = j + k = (l-|- 7 . 
Gomparing good ^ it to bad ^ 7 decay channels, here 
one defines the cooperativity parameter in a straightfor¬ 
ward way as Gatom = This is readily read as 

the cavity-to-free-space scattering ratio, since the effec¬ 
tive rate at which an excited atom emits an excitation 
into the cavity is given by it ^ For Gatom > !> the 

atom is then more likely to decay into the cavity mode 
rather than into another mode outside the cavity. In 
cavity QED, large cooperativity Gatom 1 has allowed 
for a number of key experimental demonstrations such 
as an enhancement of spontaneous emission [93] , photon 
blockade [94] and vacuum-induced transparency [95] . 

The Master equation given in Eq. (H2) describes a two- 
level system subject to purely dissipative dynamics. The 
dynamics can be fully described in terms of a set of three 
simple rate equations for the populations pj~ = {k\p\k) 
(A: = 0 , 1 ) and coherence pio = (l|/o| 0 ), summarized as 
P= {Pi,Po,Pio), 

^P= +K 0 0 \p, (H3) 

\ 0 0 -7eff/2 J 

where 7 eff = (K-fTbeph)- This allows for a simple ana¬ 
lytical solution: For example, for a system initially in the 
excited state p (t = 0 ) = | 1 ) ( 1 | it readspi (t) = exp [—kt], 
Po = 1 — exp [—Kt] and pig (t) = 0. 

Here, we aim for a theoretical description that sin¬ 
gles out the desired trajectories, where phonon emission 
through the mirrors happens first, from all others. To do 
so, we rewrite Eq.(H2) as 

P = -iHp + ipH^ + JgdP + JhdP, (H4) 

where we have defined an effective (non-hermitian) 
Hamiltonian H and jump operators according to 

H = -^7eff|l)(l| =-^(« +Edeph) |1)(1|,(H5) 

JgdP = kgdS-pS+, (H6) 

JbdP = Kbd5-p^++Edeph |l)(l|p|l)(l|. (H7) 
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Here, we have decomposed the effective decay rate k as 

k = kgd + kbd = — Kgd H-^Kbd- (H8) 

Formally solving Eq.(H4) gives 

p{t) =U it) p{Q) + [ drU {t - t) Jp{t) , (H9) 

Jo 

with the total jump operator Jp = JgoP + JoaP and 

U[t)p = (HIO) 

The exact solution given in Eq. (H9) can be iterated, giv¬ 
ing an illustrative expansion in terms of the jumps J. It 
reads 


p{t) = U(t)p(^)+ f driU {t - Ti) JU (ti) p (0) 

Jo 

rt f‘'T2 

+ / dT 2 / driU {t — T 2 ) JU {t2 — Ti) X 

Jo Jo 

JU{t^)p{0) +... 

Here, the n-th order term comprises n jumps with 
free evolution lA between the jumps. Now, we can single 
out the desired events where the first quantum jump is 
governed by jTgd- This leads to the definition 

p{t)=U {t) p (0) -f pgd (t) + Phd it) , (Hll) 

where pgd it) subsumes all desired trajectories 

/3gd(t) = [ dTiU {t - Ti) JgdU in) p (0) (H 12 ) 

Jo 

rt nT2 

-\- / dT 2 / dnU (t — T 2 ) JU (r 2 — n) X 

Jo Jo 

JgdU (ti) P (0) -I- ... 



Figure 14: (color online). Success probability pauc for a qubit 
excitation to leak through the mirror, as a function of the 
cooperativity C for e = ACbd/^gd = 0 (black solid) and e = 5% 
(blue dashed). 


We focus on a qubit, initially in the excited state, 
i.e., p(0) = |1) (1|. Using the relations U{n)pi0) = 
e- 7 e«n ^ p ( 0 ) ^ Kgde-T'=«"i | 0 ) ( 0 | and 

JgdZ^(T 2 -Ti) Jgdf^(ri)p(0) = 0, (H13) 

JbdUiT2-n)JsdUin)piO) = 0, (H14) 


the qubit’s density matrix evaluates (to all orders in J) 
to 


pit) = e '>''“*|l)(l|-bpgd(t)-bpbd(t), (H15) 
pgd(t) = M(i_e-W)|o)(o|. (H16) 

7eff 

In the long-time limit t —>■ 00 , the system reaches the 
steady state, p(t—l-oo) = pgd + Pbd) where Pgd = 
^ |0) (0|. The associated success probability Psuc = 
Tr [pgd] for faithful decay through the mirrors is then 


Psuc — 


_^gd_ 

^gd ^bd “t“ Tdeph 


1 

K _[_ 1 ^ 

^gd 4 


(H17) 


which is a simple branching ratio comparing the strength 
of the desired decay channel ~ Itgd to the undesired ones 
~ khd + Tdeph- In the limit where Kbd = 0 , i.e., k = Kgd, 
the expression for psuc simplifies to 


I /TT1 0\ 

Psnc — .. Id 
^ + 4 C 

with the usual definition found in the literature, 
C = S^/C^rdeph) = 9 '^T 2 /k; here, k = uJc/QeS = 
(hth + ljwc/Q is understood to account for thermal 
occupation of the environment hth in terms of a de¬ 
creased mechanical quality factor (compare for example 
Refs. [10, 60]). 

It is instructive to rewrite the general expression for 
Psuc given in Eq.(H17) as 


Psuc — 


1 

(1 -be) [1 -b jU] 


(H19) 


with e = Kbd/^gd- Based on this definition, it is evi¬ 
dent that two conditions need to be satisfied in order to 
reach psuc —>■ 1 in the regime where Kbd > 0 : (i) a low 
undesired loss rate, e = Kbd/^gd 1 , and (ii) high co- 

2 2 rri 

operativity, C = --Fr— = ^ 1, with n = Kgd + Kbd- 

deph ^ ® 

For an illustration, compare Fig. 14. This shows, that the 
usual definition and interpretation of the cooperativity C 
holds, provided that £ <C 1 is fulfilled. In order to quan¬ 
tify the cooperativity C for SAW cavity modes both in 
the ~ MHz (hth 1) and the ~ GHz (hth ^ 1) regime, 
in the main text we take a (conservative) estimate as 
C = ( 7 ^T 2 ( 5 /[wc (hth + 1)]. For artificial atoms (quan¬ 
tum dots, superconducting qubits, NV-centers, ...) with 
resonant frequencies GHz, at cryostatic temperatures 
this definition reduces to G « g^T 2 /K as discussed above, 
whereas for a trapped ion with u}tl2TT « ~ MHz 
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it correctly gives C « g^T^Qj {ojcfith) with a decreased 
effective quality factor Qeff = Qlnt\i [10, 96, 97]. 

Fidelity estimate. —For small errors, the expression 
given in Eq.(H19) can be approximated as Psuc ~ 1 — 
e —17(4(7). Since the absorption process is just the time- 
reversed copy of the emission process in the state transfer 
protocol for two nodes, we can estimate the state transfer 
fidelity as 7^ > Psuc x Psuc- For small infidelities, we then 
find 



where the individual errors arise from intrinsic phonon 
losses ~ £ and qubit dephasing ~ (7“^ ~ respec¬ 

tively. This simple analytical estimate agrees well with 
numerical results presented in Ref. [61], where (except 
for noise sources that are irrelevant for our problem) 
« 1 — (2/3)£/ (1 -|-£) — CC~^ with a numerical co¬ 
efficient C = O {!) depending on the specific pulse se¬ 
quence. Since e <?; 1, this relation can be simplified to 
« 1 — (2/3)£ — CC~^. For the state transfer of the co¬ 
herent superposition \%ji} = (| 0 ) — | 1 )) /v^ as described 
in detail in Appendix J, we have explicitly verified the 
linear scaling ~ e and find numerically ~ l/ 2 £ for intrin¬ 
sic phonon losses [compare also Fig. 4 in the main text, 
where T ~ 95% for e « 10% and CThuc = 0] and take 
^ £ as a simple estimate in Eq.( 8 ). Using experimentally 
achievable parameters e « 5% and C « 30, we can then 
estimate if « 90%. 


Appendix I: Effects due to the Structure Defining 
the Quantum Dots 

In our analysis of charge and spin qubits defined in 
quantum dots, we have neglected any potential effects 
arising due to the structure defining the quantum dots, 
that is (i) the heterostructure for the 2DEG and (ii) the 
metallic top gates for confinement of single electrons. In 
this Appendix, we give several arguments corroborating 
this approximate treatment, showing that the QD struc¬ 
ture does not negatively influence the cavity nor the cou¬ 
pling between qubit and cavity. 

Heterostructure. —Following the arguments given in 
Ref. [25], the 2DEG is taken to be a thin conducting 
layer a distance d away from the surface of a homoge¬ 
neous AUGai_a;As crystal with typically x « 30%. This 
treatment is approximately correct since the relevant ma¬ 
terial properties (elastic constants, densities and dielec¬ 
tric constants) of ARGai-ajAs and GaAs are very similar 
[25] . The mode-functions and speed of sound are largely 
defined by the elastic constants [25] which are roughly 
the same for both AUGai_a;As and pure GaAs; for ex¬ 
ample, the speed of the Rayleigh SAW for Alo.aGao.yAs 
is Vs « 3010m/s which differs from that of pure GaAs 
by only 5% [25] . Moreover, the numerical values for 
the material-dependent parameter H entering the ampli¬ 
tude of the mechanical zero-point motion C/q according 


to Eq.(G9) differ by 2% only [25] ; accordingly, the es¬ 
timate for our key figure of merit C/q should be rather 
accurate. Also, the piezo-electric coupling constants are 
rather similar, with ei 4 « 0.15C/m^ for pure GaAs and 
ei 4 « 0.145C/m^ for Alo. 3 Gao. 7 As [24, 25] yielding an 
accurate estimate for ^o and ^oi respectively. Lastly, 
the heterostructure is not expected to severely affect the 
(5-factor, since very high (5-values reaching Q > 10^ 
have been observed in previous SAW experiments on 
AlN/diamond heterostructures [53, 54], where the differ¬ 
ences in material properties are considerably larger than 
for the heterostructure making up the 2DEG. 

Top gates. —For the following reasons, we have disre¬ 
garded effects due to the presence of the metallic top 
gates: (i) In Ref. [27], a closed form analytic solution for 
the piezoelectric potential (/) (x, t) accompanying a SAW 
on the surface of a GaAs/AUGai_a;As in the presence 
of a narrow metal gate has been obtained. In partic¬ 
ular, it is shown that (j) (x, t) is screened right below 
the gate, but remains practically unchanged with respect 
to the ungated case outside of the edges of the gate. 
Since the QD electrons are confined to regions outside 
of the metallic gates, they experience the piezoelectric 
potential as calculated for the ungated case (see Ap¬ 
pendix B for details); therefore, our estimates—where 
this screening effect has been neglected—remain approx¬ 
imately valid, (ii) In Ref. [46], the coupling of a traveling 
SAW to electrons confined in a DQD has been experi¬ 
mentally studied. Here, in very good agreement to the 
experimental results, the potential felt by the QD elec¬ 
trons has been taken as Vpe sin (kl) (where I is the 
lithographic distance between the dots) confirming the 
sine-like mode profile as used in our estimates. More¬ 
over, with I = 220nm and A = l.d/rm the parameters 
used in this experiment perfectly match the ones used 
in our estimates. Intuitively (in the spirit of the stan¬ 
dard electric dipole approximation in quantum optics), 
since A ^ Z, the SAW mode cannot resolve the dot struc¬ 
ture and thus remains largely unaffected, (iii) In Ref.[16[, 
single phonon SAW pulses have been detected via a sin¬ 
gle electron transistor (SET) directly deposited on the 
GaAs substrate with time-resolved measurements clearly 
identifying the coupling as piezo-electric. Similarly to a 
QD, the SET is defined by metallic gates. Here, a re¬ 
lation between vertical surface displacement and surface 
charge induced on the SET is theoretically derived. Us¬ 
ing standard tabulated parameter values for GaAs and 
neglecting any effects due to the presence of the metallic 
gates, very good agreement with the experimental results 
is achieved. In particular, based on the results given in 
Ref. [16], a straightforward estimate gives Uq « 30am for 
the rather large cavity with A « 10®^m^, whereas Eq.(4) 
yields a smaller, conservative estimate Uq « 2am, due 
to the averaging over the quantization area A. (iv) The 
( 5 -factor of the cavity is not expected to be severely af¬ 
fected by the presence of the metallic gates since metal¬ 
lic A1 cladding layers have been used on a GaAs sub¬ 
strate to show basically dissipation-free SAW propaga- 
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tion over millimeter distances [38, 39]. (v) Finally, there 
is a large body of previous theoretical works on electron- 
phonon coupling in gate-defined QDs (see for example 
Refs. [24, 87, 88, 98[) where any effects due to the struc¬ 
ture defining the QDs have been neglected as well. As 
a matter of fact, our description for the electron-phonon 
coupling emerges directly from these previous treatments 
in the limit where the continuum of phonon modes is re¬ 
placed by a single relevant SAW cavity mode (similar to 
cavity QED, other bulk modes are still present and con¬ 
tribute to the decoherence of the qubit on a timescale 
T 2 ). For example, the piezoelectric electron-phonon in¬ 
teraction is given in Ref.[98[ as 


T^GaAs 


/^E 


2 pVv.k 


M, 


0'k,fi + fl 


t 

— k,ii ’ 


(II) 


where a\ creates an acoustic phonon with wave vector 
k and polarization p and Aik refers to the matrix ele¬ 
ment for electron-phonon coupling; for free bulk modes, 
Aik = Ha (*l \j) 4adja- In agreement with 

our notation, the coupling constant is /? = eei 4 /e and 
the square-root factor can be identified with Uq. Re¬ 
placing the sum by a single-relevant cavity mode a, we 
recover the Hamiltonian describing the cavity-qubit cou¬ 
pling with g ~ pUo = e4>Q\ compare Eq.(Fl). 


Appendix J: State Transfer Protocol 

In this Appendix, we provide further details on the 
numerical simulation of the state transfer protocol as de¬ 
scribed by Eqs.(9) and (10), in the presence of Marko¬ 
vian noise. In line with previous theoretical studies [97[, 
we show that the simple approximate Markovian noise 
treatment results in a pessimistic estimate for the noise 
transfer fidelity A. 

We first provide results of the full time-dependent nu¬ 
merical simulation of the cascaded Master equation given 
in Eqs.(9) and (10), including an exponential loss of co¬ 
herence for Fdeph > 0; see Fig. 15. In contrast to the 
non-Markovian noise model discussed in the main text, 
the qubits are assumed to be on resonance throughout 
the evolution, that is <5^ = 0 (i = 1,2), but experience 
undesired noise as described by 

-^noise/^ — ^ ^ [^3-2] p ~\~ Tfieph (Ji) 

i=l,2 i-1,2 

Here, the second term refers to a standard Markovian 


pure dephasing term that leads to an exponential loss 
of coherence ~ exp (—Fdepht/S). As a reference we 
also show the results for the ideal, noise-free scenario 
(£noise = 0), where perfect state transfer is achieved [20[. 
The results of this type of time-dependent numerical sim¬ 
ulations are then summarized in Fig. 16. For optimized, 
but experimentally achievable parameters T 2 « 3^s 
[68[, and accordingly FdeDh/Ksd ~ 3%, we then obtain 



time Kgdt 


Figure 15: (color online). Numerical simulation of state 
transfer for two different initial states |l)j (red solid line) 
and (|0)j^ — |1)]^) /y/2 (blue dashed line) in the presence of 
Markovian noise. Black (dash-dotted) curves refer to the 
ideal, noise-free scenario (Tnoise = 0) where perfect transfer is 
achieved, while colored curves take into account decoherence 
processes, (a) Transfer fidelity A; (b) Excited-state occupa¬ 
tion for first and second qubit. Numerical parame¬ 

ters. gj_ (t ^ 0 ) — Atgd, ^bd/^gd — 5?^ and F^eph/^gd — 59o. 


A « 0.85, in the presence of realistic undesired phonon 
losses Khd/ngd = 5%. 

This shows that our non-Markovian noise model yields 
even higher state transfer fidelities than the Markovian 
noise model. Intuitively, this can be readily under¬ 
stood as follows: The simple Markovian noise model 
gives a coherence decay ~ exp [—t/T 2 ), whereas our non- 
Markovian noise model yields ~ exp(—There¬ 
fore, for Markovian noise one can estimate the dephasing 
induced error on the relevant timescale for state trans¬ 
fer ~ as ^ K~^/T 2 ~ Fdeph/^ ^ 1) whereas non- 
Markovian noise leads to a considerably smaller error 
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